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In the initial years of operation, experiments at the Relativistic Heavy Ion Collider 
\ (RHIC) have identified a new form of matter formed in nuclei-nuclei collisions at 

energy densities more than 100 times that of a cold atomic nucleus. Measurements 
and comparison with relativistic hydrodynamic models indicate that the matter 
thermalizes in an unexpectedly short time, has an energy density at least 15 times 
larger than needed for color deconfinement, has a temperature about twice the crit- 
ical temperature predicted by lattice QCD, and appears to exhibit collective motion 
with ideal hydrodynamic properties - a "perfect liquid" that appears to flow with a 
near-zero viscosity to entropy ratio - lower than any previously observed fluid and 
perhaps close to a universal lower bound. However, a fundamental understanding 
of the medium seen in heavy- ion collisions at RHIC does not yet exist. The most 
important scientific challenge for the field in the next decade is the quantitative 
exploration of the new state of nuclear matter. That will require new data that 
will, in turn, require enhanced capabilities of the RHIC detectors and accelerator. 
In this report we discuss the scientific opportunities for an upgraded RHIC facility 
- RHIC II - in conjunction with improved capabilities of the two large RHIC detec- 
tors, PHENIX and STAR. We focus solely on heavy flavor probes. Their production 
rates are calculable using the well-established techniques of perturbative QCD and 
their sizable interactions with the hot QCD medium provide unique and sensitive 
measurements of its crucial properties making them one of the key diagnostic tools 
available to us. 
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1 Introduction 

During the last 7 years, heavy-ion experiments at the Relativistic Heavy-Ion 
Collider (RHIC) have recorded a wealth of data in Au+Au, d+Au, Cu+Cu, 
and pp collisions at a variety of energies, ranging from ^snn = 19.6 GeV 
to the highest available Au+Au energy, 200 GeV. It is at these high energies 
that QCD predictions of new phenomena come into play under conditions 
where, over nuclear volumes, the relevant degrees of freedom are expected 
to be those of quarks and gluons rather than of hadrons, the realm of the 
quark-gluon plasma. 

Measurements from the four RHIC experiments, BRAHMS, PHENIX, PHO- 
BOS and STAR, have revealed compelling evidence for the existence of a new 
form of nuclear matter at extremely high densities and temperatures [1-4]. 
Detailed analyses of these data also make it clear that this hot, dense medium 
has surprising properties. 

The properties of the medium are those of a strongly coupled plasma, or 
sQGP, that behaves like a "perfect liquid" flowing with a near-zero viscosity 
to entropy ratio [5] . The RHIC observations have spurred significant advances 
in theory. However, a fundamental understanding of the medium seen in heavy- 
ion collisions at RHIC does not yet exist. It requires new data that, in turn, 
necessitate enhanced capabilities of the RHIC detectors and accelerator. A 
detailed plan is being developed by BNL to implement these upgrades in 
collaboration with the RHIC scientific community. 

The main focus of this report is to outline the scientific opportunities in 
the heavy flavor sector provided by upgrades of the two large RHIC detec- 
tors 1 , PHENIX and STAR, in conjunction with an upgrade of the accelera- 
tor/collider facility, referred to as RHIC II. The detector upgrades will improve 
the acceptance, particle identification and secondary vertex detection capabil- 
ities of PHENIX and STAR. The RHIC II accelerator luminosity upgrade was 
originally proposed to be done using electron cooling. It is now planned to 
use a scheme that employs stochastic cooling and other upgrades to achieve 
for Au+Au a factor of 5 increase in luminosity over present capabilities. This 
will provide about 70% of the average luminosity expected from electron cool- 
ing, but much sooner and at a much lower cost. The accelerator upgrades will 

1 The two smaller experiments BRAHMS and PHOBOS were decommissioned in 
2006. 
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include a new ion injector, EBIS, which will provide high-intensity beams of 
nuclei as massive as uranium. 

This report is the result of the collaboration and research efforts of a RHIC- 
wide Heavy Flavor Working Group. It provides a comprehensive overview of 
the physics questions than can be addressed by studies of open charm, open 
bottom and quarkonia at RHIC II. It also includes a detailed assessment of the 
accelerator and detector capabilities required to carry out these measurements 
with sufficient precision to resolve many of the outstanding issues by providing 
detailed results with which to make thorough comparisons to current and 
future theoretical calculations. 

This report is organized as follows. After a general introduction to heavy 
flavor physics, section 2 discusses the detector upgrade program at RHIC. The 
projected yields of various heavy flavor measurements that can be achieved 
utilizing these upgrades and the higher RHIC II luminosities are discussed 
in section 3. In sections 4 and 5 we present a more detailed discussion of 
the motivation for studying open heavy flavor and quarkonia in heavy-ion 
collisions, respectively. We also include a review of the current theoretical and 
experimental status as well as the proposed experimental program. In section 
6 we review the relationship between heavy flavor physics at RHIC II and the 
LHC. We conclude in section 7. 



1.1 Motivation 

Because charm and bottom quarks are massive, they are produced almost 
exclusively in the initial parton-parton interactions in heavy-ion collisions at 
RHIC energies. In the absence of any nuclear effects, the heavy flavor cross 
sections in A + A collisions at RHIC would simply scale with the number of 
binary collisions. Thus departures from binary scaling for heavy flavor produc- 
tion in A + A collisions provide information about nuclear effects. These can 
be divided into two categories: effects due to embedding the colliding partons 
in a nucleus (cold matter effects) and effects due to the large energy density 
in the final state. The main focus of the heavy flavor program at RHIC is to 
investigate the properties of the dense matter produced in A + A collisions by 
studying its effects on open heavy flavor and quarkonium production. This in 
turn requires a detailed understanding of cold matter effects so that they can 
be unfolded from the dense matter effects. 

The program thus requires detailed measurements and calculations of pp and 
p + A heavy flavor cross sections to characterize the cold matter effects, if 
we are to quantify the differences between QGP and non-QGP effects. Up-to- 
date benchmark calculations of the total open heavy flavor (charm and bottom 
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hadrons) and quarkonium (J/ip and T families) yields and spectra are imper- 
ative. Cold matter effects that need to be included are nuclear shadowing, 
for both open heavy flavor and quarkonium production, and nuclear absorp- 
tion of quarkonium. Recent calculations of charm and bottom production to 
Fixed-Order Next-to-Leading Logarithm (FONLL) in pp collisions have been 
published [6], along with a discussion of the inherent theoretical uncertainties 
[7] and reference calculations of heavy quark, heavy flavor meson and their 
decay lepton spectra [6]. Similar calculations have been made for quarkonium 
production, including studies of shadowing and absorption effects as a function 
of rapidity and centrality in d+Au [8] and A + A [9] collisions at RHIC. 

A number of dense matter effects on heavy flavor production have been pre- 
dicted. Some of these do not change the total cross section but, instead, modify 
the pr spectra of heavy flavor hadrons and their decay products. Heavy quark 
energy loss [10-14] by collisional and radiative processes steepens the p T dis- 
tribution relative to that in pp collisions. On the other hand, random pt kicks 
result in transverse momentum broadening, increasing the average pr in both 
cold nuclear matter [15] and in passage through hadron bubbles in the mixed 
phase of a QGP [16]. If the medium surrounding the heavy quarks after pro- 
duction exhibits collective motion, such as transverse flow [17, 18], the low p T 
heavy quarks (pr < m) may be caught in this flow. Strong effects of energy 
loss [19,20] on heavy flavor decays to electrons and charm flow [20] have al- 
ready been seen in Au+Au collisions at RHIC. Studying heavy flavor energy 
loss using single electrons requires being able to separate electrons from c and 
b decays since the large bottom and charm quark mass difference suggests 
that bottom quark energy loss is weaker than that for charm [10]. Some QGP 
studies require accurate baseline determinations of the total heavy flavor cross 
sections to interpret other effects. For example, if more than one cc pair is pro- 
duced in an A + A event, uncorrelated c and c quarks might coalesce to form 
a J/^i in a QGP [22-25]. The total cc yield is needed to normalize the J/ip 
production rate from this process. 

Suppression of J/ip production was one of the most exciting proposed QGP 
signatures at the CERN SPS [26]. This J/ip suppression was predicted to oc- 
cur due to the shielding of the cc binding potential by color screening, leading 
to the breakup of the quarkonium states, first the Xc and ip', and finally the 
J I ip itself as the temperature increases [27, 28] . The QGP suppression may 
not be so simple, as lattice gauge theory studies of the J/ip spectral function 
above the critical temperature for deconfinement, T c , attest. The J/ip may 
exist as a bound state for temperatures considerably larger than T c [29]. How- 
ever, the J ftp may instead be dissociated by hot thermal gluons in medium 
[30] before it could be suppressed by color screening. Secondary quarkonium 
production from uncorrelated QQ pairs, either in the plasma phase [23, 25, 31- 
33] or in the hadron phase [34,35], could counter the effects of suppression, 
ultimately leading to enhanced quarkonium production. Such secondary J/-0's 



5 



would have different kinematic distributions than those from the initial pro- 
duction. Because the underlying cc distribution falls rapidly with p?, the pt 
distribution produced by coalescence will be softer. If the underlying cc distri- 
bution peaks at midrapidity, the J/ip rapidity distribution from coalescence 
will be narrower than that produced in the primordial collisions. The coa- 
lescence rapidity distribution should be calculated with shadowing effects on 
the underlying cc distribution taken into account since these can cause the cc 
distribution to flatten in more central A + A collisions [9]. Elliptic flow is also 
expected to affect quarkonium as well as open heavy flavors [17, 18]. 

With higher luminosity at RHIC, the T yields could also be measured ac- 
curately. Since the T radius is smaller than that of the J/ip [28], direct 
color screening in the QGP would not occur until much higher temperatures. 
The higher mass bottomonium states, however, would likely be suppressed 
at RHIC, as would the Xc and ip' in the charmonium family. The feed down 
structure is more complicated for the T since there are three S states (T, T' 
and T") and two sets of P states (xbi and Xb2) below the BB threshold. The 
T family suppression should be measurable over a large px range, with QGP 
suppression possible on the T' and T" up to p T ~ 40 GeV/c at the LHC [36]. 
Because of the small number of bb pairs produced at RHIC, bottomonium 
formation by coalescence of unrelated pairs should be negligible. 



1.2 Overview of results from the heavy flavor program at RHIC 

Heavy flavor measurements capable of discriminating between theoretical mod- 
els need large integrated luminosity. In RHIC runs so far, useful data sets have 
been acquired at 200 GeV for pp, d+Au, Cu+Cu and Au+Au collisions. The 
Run 6 pp data are not fully analyzed yet, but a sizable data set is available 
from Run 5. Recent runs for two species greatly increase the statistical reach 
over earlier runs: Au+Au (Run 7) and d+Au (Run 8). The data are still being 
analyzed, but some preliminary results are available from Run 7. 

The pp data collected to date provide an essential reference for the heavy- 
ion program in the form of the underlying heavy flavor production rates as 
functions of rapidity and pt- Equally essential, the data from d+Au collisions 
provide baseline information about cold nuclear matter effects which must also 
contribute to heavy flavor production in heavy-ion collisions. The analyzed 
d+Au data from Run 3 have limited statistical precision, but they provide 
useful tests of models that include the effects of shadowing on heavy flavor 
production and of J/tp absorption in cold nuclear matter [8]. 

The easiest way by far to measure open heavy flavor yields in heavy-ion col- 
lisions at RHIC is via the semileptonic decays of D and B mesons. Two 
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very striking and unexpected results have already been seen by studying de- 
cay electrons from open heavy flavor in heavy-ion collisions at RHIC. The 
first is the observation that the nuclear modification factor, Raa(v,Pt) — 
(daAA/dprdy) / ({TAB)d°'pp/ dprdy) , for electrons from open heavy flavor de- 
cays shows very strong suppression in central Au+Au collisions [19,20], sim- 
ilar to that seen for pions. The second striking result is that the elliptic flow 
parameter, V2, of electrons from open heavy flavor decays appears to favor 
charm quark flow at low p T [20]. Until recently, it had been expected that 
heavy quark energy loss would be considerably smaller than that for light 
quarks due to interference effects [10]. Generating the necessary energy loss 
for charm and bottom quarks with realistic gluon densities in the material is 
a major challenge for models [10,37]. The relatively large t>2 values at low pr 
imply at least some degree of charm quark equilibration with the medium. 
This also implies very strong interactions of charm quarks with the medium 
at lower p T [17, 18]. 

A serious shortcoming of open heavy flavor measurements employing semilep- 
tonic decays is the difficulty of separating the contributions to the lepton 
spectra from charm and bottom decays. Perturbative calculations of the rel- 
ative contributions from charm and bottom as a function of Pt have large 
theoretical uncertainties: at midrapidity the crossover point at which bottom 
decays become dominant is between 3 and 9 GeV/c. Recently, however, the 
ratio of charm to bottom contributions has been extracted by STAR and 
PHENIX from pp data, although with still large uncertainties. The STAR 
measurements [38] cover an electron p? of 3 — 9 GeV/c and are based on small 
azimuthal angle correlations between the decay electrons and hadrons as well 
as correlation between electrons and identified D° mesons. The PHENIX bot- 
tom/charm ratio as a function of electron p? [39] was inferred from the shape 
of the electron-hadron correlation function. The STAR and PHENIX ratios 
are in good agreement and indicate that bottom becomes the dominant con- 
tribution to the electron spectra at pr ~ 3.5 GeV/c (see Fig. 6). A separate 
PHENIX measurement of the bottom and charm cross sections was also ob- 
tained from fits to the dielectron invariant mass spectrum, after removal of 
all other contributions [40]. The results are in agreement with the cross sec- 
tions inferred by PHENIX by combining the correlation measurements with 
electron pt distributions. 

PHENIX has also measured open heavy flavor yields at forward rapidity, 
y = 1.65, using single muons [41]. The charm cross section is da/dy = 
0.145 ± l.l%+4gg% lib. This is consistent with the measured values at y = 
but the precision is currently insufficient to meaningfully define the rapid- 
ity dependence. It is likely that the displaced vertex measurement from the 
planned vertex detector upgrades will be needed to determine the rapidity 
distribution. 



7 



The first high statistics charmonium results for heavy-ion collisions at RHIC 
are now available [42,43], along with high statistics pp reference data [44]. 
These include Au+Au and Cu+Cu results for the J/ip Raa as a function of 
the number of participant nucleons, N pait , in the rapidity intervals \y\ < 0.35 
and 1.2 < \y\ < 2.2. A striking feature of the Au+Au J/ip data is that, for 
-Wpart > 150, the suppression is considerably stronger at forward rapidity than 
at midrapidity. Comparison with existing models at midrapidity shows that 
cold nuclear matter baseline calculations [9] which approximately reproduce 
the RHIC d+Au J/ip rapidity distributions [8] somewhat underpredict the 
suppression observed in Cu+Cu and Au+Au collisions. On the other hand, 
several suppression models [25,45,46] which were successful in describing J/ip 
suppression at the SPS are found to strongly overpredict the suppression at 
RHIC. Models which incorporate strong suppression combined with J/ip coa- 
lescence from uncorrelated cc pairs seem to agree best with the data although 
the existing models slightly underpredict the suppression. A major source of 
uncertainty in all of this is due to the poor precision of the existing d+Au 
data, leading to poor constraints on the baseline cold nuclear matter effects. 

A recent reanalysis of the Run 3 d+Au J/ip data [47] used the higher yield 
Run 5 pp data as the reference for calculating Raa- Theory calculations have 
been fit to the reanalyzed data to explore how well the data constrain cold 
nuclear matter effects in heavy-ion collisions. They find that the constraints 
from the Run 3 d+Au data are insufficient for firm conclusions to be drawn 
about additional hot matter effects in central heavy-ion collisions [43,47]. The 
newly obtained Run 8 d+Au data set, approximately 20 times the yield of 
Run 3, will provide much improved statistical precision. 

In the last few years, theorists have begun exploring the consequences of J / ip 
coalescence on observables other than the centrality dependence of the nuclear 
modification factor [22,48,49]. This work has led to the prediction that J/ip's 
formed by coalescence of uncorrelated cc pairs will have narrower rapidity and 
Pt distributions due to the presumed shape of the underlying charm quark dis- 
tributions. The coalescence contribution to J/ip production will cause many 
observables to change with centrality, including the rapidity and pt depen- 
dence of Raa, the shape of the pt distribution (quantified by the average ptj,, 
(Pt))> an d the J/ip elliptic flow parameter, t> 2 . Predictions of (p^) as a function 
of centrality have been made for Au+Au and Cu+Cu collisions with and with- 
out coalescence [22,48,49]. These predictions, compared to the Au+Au and 
preliminary Cu+Cu data, favor calculations including coalescence. The coales- 
cence contributions predicted for central Au+Au (as well as central Cu+Cu) 
collisions [22] are also qualitatively expected to narrow the J/ip rapidity dis- 
tributions if the underlying charm distributions are peaked at midrapidity. 
The Au+Au data [42] show some narrowing of the rapidity distribution for 
the most central collisions although the reduction in RMS is only about 2a. 
Work is still needed to quantify both the theoretical predictions and the ex- 
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perimental observables. A first measurement of the J/ip v 2 may come out of 
the Run 7 Au+Au data, albeit with large uncertainties. 

The existing, but not yet analyzed Au+Au data from Run 7 will quantita- 
tively improve the measurements of many heavy flavor observables. The J/ip 
Raa versus centrality and (p\) as well as Raa and v 2 measurements of charm 
and bottom semileptonic decays to single electrons will all improve signifi- 
cantly, allowing more stringent model tests. Measurements of other observables 
will be qualitatively improved. Examples are: definitive v 2 measurements from 
semileptonic decays at intermediate to high p T where we might hope to see 
the transition from charm to bottom dominance and flow to non-flow; a pos- 
sible first J /ip v 2 measurement; substantially improved measurements of J/ip 
Raa{u) to quantify the coalescence contribution; and improved measurements 
of J /if) Raa{pt), invaluable for understanding coalescence and formation time 
effects. 

However it is clear that the RHIC heavy flavor program is now limited by 
the capabilities of the accelerator and the detectors. The accelerator upgrades 
planned over the period 2009-2013 to produce a factor of 5 greater luminosity 
at RHIC II, combined with the detector upgrades in place by that time, will 
be required for the heavy flavor program at RHIC to move to the next level, 
as described below. 

1.3 Overview of the proposed heavy flavor program at RHIC II 

The increase in luminosity, combined with the increased capabilities of the 
upgraded PHENIX and STAR detectors, will make it possible to add many 
important new probes to the heavy flavor program at RHIC. 

One of the most powerful benefits of the luminosity upgrade will be the ability 
to measure yields of the excited charmonium states: the ip' and Xc- Lattice 
calculations predict much smaller melting temperatures for the ip' and Xc than 
for the more tightly bound J /if). Thus these excited states should not be able 
to exist in the QGP at RHIC and comparison of the tp' and Xc yields to the 
J /ip yield as a function of centrality would be a direct test of deconfinement. 

Testing models in which the observed J/ip yield in heavy-ion collisions is due 
to competition between gluon dissociation and coalescence formation in the 
QGP requires very high luminosity. Tests of charm coalescence models include 
measuring J/ip v 2 as a function of px, J/tp Raa to much higher p T to follow the 
trends of suppression as the J/ip formation time approaches the QGP crossing 
time, and J/ip polarization as a function of collision centrality. The rapidity 
and pt dependence of Raa as functions of ^s^n and centrality, requiring 
sufficient luminosity for precision measurements at multiple energies, is not 
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possible on a reasonable time scale at the present RHIC luminosity. 

The detailed study of bottomonium states, the T family, is only possible at 
RHIC II luminosities. Like the charmonium states, the dissociation tempera- 
tures of the bottomonium states depend on the binding energies. There are, 
however, two important differences. First, the bottomonium binding ener- 
gies, particularly that of the T(1S), are higher so that they should dissociate 
at higher temperatures. Only the higher-lying bottomonium states are thus 
likely to break up at RHIC energies. Second, the bb production rate in cen- 
tral Au+Au collisions is only ~ 0.05 pairs per collision, making coalescence 
production of bottomonium much less likely. Thus bottomonium production 
at RHIC II will provide a very different window on color screening effects 
than charmonium production. The bottomonium yields at RHIC II should be 
sufficient for measurements of Raa as a function of centrality in heavy-ion 
collisions for the three T S states. The T yields at RHIC II and at the LHC 
will not be sufficient for t> 2 or polarization measurements. 

As mentioned earlier, measurements of semileptonic open heavy flavor decays 
at RHIC have already produced strikingly different results than expected. 
The strong suppression in Raa coupled with the large v 2 suggest very large 
heavy quark energy loss in the medium. The fact that these semileptonic 
decay spectra contain both charm and bottom contributions remains a sig- 
nificant complication, in spite of recent successes by STAR and PHENIX at 
measuring the ratio of b to c quark contributions in the semileptonic decay 
spectrum of pp collisions. The separation of open charm and bottom can be 
done in several ways. Charm can be observed via hadronic D° — > K ± n T and 
— > K ± 7t ± tt t decays, as STAR does. Precise Raa and i>2 measurements are 
difficult in this channel. Since these events cannot be triggered, they must be 
extracted from a minimum bias data set that samples only a small fraction 
of the available luminosity. The combinatorial background is also very large, 
making statistical precision difficult. The addition of a displaced vertex mea- 
surement in STAR will dramatically reduce the combinatorial background but 
there is still no trigger for these decays. At RHIC II luminosity, bottom can 
be observed very cleanly in both PHENIX and STAR via B — > J / ip X decays 
using displaced vertices, providing good measurements of the bb cross section 
and bottom quark Raa- However those yields will certainly be too small for 
v 2 measurements at RHIC II or the LHC. Finally, the combination of RHIC 
II luminosity with a displaced vertex measurement should allow statistical 
separation of the charm and bottom contributions to the semileptonic decay 
spectra, taking advantage of the different c and b quark decay lengths. Such 
semileptonic decay measurements, while less clean than the direct D and B 
decay measurements, have the advantage of much larger yields so that separate 
charm and bottom v 2 measurements should be possible. 

Independent measurements of open charm and bottom Raa and vi to high pt 
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will be a very important capability at RHIC II. At low p?, these measurements 
reflect the degree of heavy quark thermalization in the medium. At high p?, 
they probe the energy loss of heavy quarks in the medium, providing an inde- 
pendent measurement of the initial energy density relative to the light quark 
energy loss measurements. The thermalization and energy loss mechanisms at 
low and high px respectively may be quite different due to possible resonance 
scattering at low pr- 



1.4 Overview of the relationship of RHIC II to the LHC program 

The heavy flavor production cross sections are significantly higher at the LHC 
than at RHIC since the per nucleon Pb+Pb energy at the LHC is a factor of 
27.5 higher than the maximum per nucleon Au+Au energy at RHIC. The cc 
and bb cross sections increase by factors of 15 and 100 respectively [15] while 
the J/vp and T cross sections increase by factors of 13 and 55 respectively 
[50]. But, because of the 10 times higher average luminosity and three times 
longer heavy-ion runs, the Au+Au integrated luminosity at RHIC II will be 
much higher than for Pb+Pb at LHC. Therefore the heavy flavor yields per 
year are expected to be similar at the two facilities. 

At y^s = 200 GeV, bottom decays to leptons begin to dominate the single 
electron spectrum at pr ~ 4 GeV/c. As the collision energy increases, the 
lepton spectra from B and D decays move closer together rather than further 
apart [15]. Thus, the large increase in the bb cross section relative to cc does 
not make single leptons from B and D decays easier to separate. Preliminary 
calculations show that the B — > e decay does become larger than that of 
D — > e, but at pr > 10 GeV/c. The yields from the two lepton sources differ 
by less than a factor of two up to pr ~ 50 GeV/c in the range \y\ < 1. Thus 
interpretation of single lepton results on heavy flavors will be more difficult at 
the LHC. Other means of separating charm and bottom must be found. ALICE 
can reconstruct hadronic D° decays from p? ~ to px ~ 25 GeV/c [51] but, 
like STAR, will have to rely on minimum bias data for these measurements 
because of the lack of a trigger. While it is not yet clear what CMS and ATLAS 
will do to reconstruct charm, they should be able to make 6-jet measurements, 
similar to the Tevatron. One way that B mesons can be measured at the LHC 
is through their decays to J/ip, as discussed further below. It has also been 
suggested that the BB contribution to the dimuon continuum, the dominant 
contribution above the T mass, can be used to measure energy loss [52] . That 
channel would be fairly clean at the LHC but more difficult at RHIC. 

The RHIC II upgrades and the high LHC energies make detailed studies of T 
production and suppression possible. At the LHC, higher initial temperatures 
make T suppression more likely than at RHIC II. But the higher bb production 
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rate (~ 5 per central Pb+Pb collision) means that, unlike RHIC, significant 
coalescence contributions to T production may be expected at the LHC Thus 
measurements at the two energies complement each other. At RHIC II, it is 
likely that PHENIX will be able to measure and resolve the three T S states. 
STAR and ALICE will have similar T yields but the STAR mass resolution 
will require fitting to extract yields. The CMS detector at the LHC has suf- 
ficient mass resolution to separate all three T S states. The T states can be 
measured to pt ~ at all LHC detectors. Only ALICE will be able to mea- 
sure J/ip production to p T ~ without a special trigger [50] since CMS and 
ATLAS require high single muon p T so that typically only J / ip with p T > sev- 
eral GeV/c are accepted. (However, CMS is working on a higher- level trigger 
to measure lower pt J/ip [53].) The larger bb cross section at the LHC means 
that J /ip production from B — > J/ipX cannot be neglected. These decay J/ip 
should be separable from the initial production using displaced vertices [50] . 



2 Detector upgrade program at RHIC 

Both PHENIX and STAR have extensive upgrade programs underway that 
are extremely important for the heavy flavor program. The upgrades that are 
most relevant to heavy flavor measurements are described here. The impact 
on the heavy flavor program of these detector upgrades, in combination with 
the RHIC II luminosity increase, will be discussed in sections 4 and 5. 



2.1 PHENIX upgrades 

Several PHENIX detector upgrades that greatly enhance the heavy flavor 
capability of the experiment are expected to be available in the RHIC II 
time frame. The most important upgrades for the heavy flavor program will 
be the barrel [54] and endcap [55] Silicon Vertex Detectors, the Nose Cone 
Calorimeter [56] and the Muon Trigger Upgrade [57]. The central region of 
the PHENIX detector, after installation of the silicon trackers and the Nose 
Cone Calorimeter, is shown in Fig. 1. The pseudorapidity and azimuthal angle 
coverages of the new detectors are illustrated in Fig. 2. 

The Silicon Vertex Detector (SVTX) consists of a central barrel [54] and two 
endcap detectors [55], as shown in Fig. 1. The SVTX barrel will have a dis- 
placed vertex resolution of ~ 50 /zm while the endcap resolution is ~ 90 — 115 
fim. Together, this inner tracking system provides full azimuthal coverage 
over \r)\ < 2.4. The SVTX will tag heavy flavor decays using displaced ver- 
tices by connecting to tracks in both the central and muon arms, improving 
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Fig. 1. The central region of the PHENIX detector after the addition of the barrel 
and endcap silicon vertex detectors and the Nose Cone Calorimeter. 




Fig. 2. The pseudorapidity and azimuthal angle coverage of the PHENIX barrel and 
endcap silicon vertex detectors and the Nose Cone Calorimeter (NCC) [56]. Two 
areas of the central arms that that provide hadron identification to high pyare also 
shown. 

the quarkonium invariant mass resolution and reducing backgrounds to heavy 
flavor measurements. In the muon arms, a loose displaced vertex cut will 
eliminate most muon tracks from light hadron decays and a very tight cut, 
~ 2cr where a is the resolution of the displaced vertex measurement, will elim- 
inate most punch-through hadrons. The displaced vertex measurement will 
greatly enhance D° — > K^-n^ measurements in the central arms, presently 
very difficult in PHENIX, by reducing the contribution to the combinatorial 
background both from prompt tracks (by using a tight vertex cut) and light 
meson decay tracks (by using a loose cut of ~ 1 cm). A loose displaced ver- 
tex cut will also reduce high background tracks in the central arms due 
to misidentified light hadron decays. In addition to identifying semileptonic 
heavy flavor decays, displaced vertex measurements can help identify J/-0's 
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from B meson decays since all other J/^'s are prompt. 

The SVTX barrel is presently under construction. It will consist of four con- 
centric silicon layers. The two inner layers, at radii of 2.5 and 5.0 cm, consist of 
pixel detectors with a segmentation of 50 /im by 425 /xm. The outer two layers, 
with radii of 10 and 14 cm, consist of 80 fim by 3 cm strips. The occupancy 
of the inner layer will be about 4.5% in central Au+Au collisions. The SVTX 
barrel produces a dramatic improvement in high track resolution in the cen- 
tral arms. The PHENIX Drift Chamber is outside the magnetic field so that, 
in the present momentum measurement, there is no information about the 
initial azimuthal angle, 0, of the track. The momentum is calculated from the 
difference between the angle of the track after passing through the magnetic 
field and that from the vertex position to the Drift Chamber. This difference 
is only ~ 40% of the total deflection. By adding a precise measurement of 
the initial direction, the SVTX barrel measures the full deflection directly, 
decreasing the momentum resolution by a factor of ~ 2.5, greatly improving 
the T invariant mass resolution. Installation of the barrel is expected starting 
in 2009. 

The forward silicon detector endcaps will consist of four silicon mini-strip 
planes. The mini-strips have 75 /im pitch in the radial direction and lengths 
in the <fi direction varying from 2.8 mm to 12.1 mm, depending on the polar 
angle. The maximum occupancy per strip is estimated to be less than 2.8% 
in central Au+Au collisions. The displaced vertex resolution of 90 — 115 /im, 
depending on the number of layers of silicon traversed by the track, should 
be compared to a mean vertex displacement of 785 /im for the boosted open 
charm muons. A prototype covering about 1/8 of one muon arm is presently 
under construction. 

The PHENIX Nose Cone Calorimeters (NCCs) [56] , tungsten-silicon calorime- 
ters, will replace the two central arm magnet nose cones, and will cover 
0.9 < \r]\ < 3.5. The simulated energy resolution for photons is ~ 27%/\/E 
GeV. The Nose Cone Calorimeters will contain both electromagnetic and 
hadronic calorimeter sections. The electromagnetic calorimeter will contain 
a pre-shower detector and a shower-max detector designed to discriminate be- 
tween individual electromagnetic showers and overlapping photons from high 
momentum 7r° decays. The pre-shower and shower-max detectors are expected 
to resolve showers with separations down to 2 and 4 mm, respectively. The 
NCCs should thus have good acceptance for Xc J/i> + 7 decays with the 
J/ip detected in the muon arms. 

The muon trigger upgrade [57] is required for PHENIX to be able to take 
complete advantage of the RHIC II luminosity upgrade for muon arm mea- 
surements. The current muon arm level- 1 heavy vector meson triggers have 
sufficient rejection capability to handle Au+Au collision rates of up to ~ 20 
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kHz and pp collision rates of up to ~ 0.5 MHz. The muon trigger upgrade 
adds three layers of Resistive Plate Chamber (RPC) detectors, with two di- 
mensional (6, 0) readout, in each muon arm. These layers follow the design 
of the CMS muon trigger at the LHC but the cathode pad segmentation is 
optimized for PHENIX. The front end electronics and trigger processors will 
be developed within PHENIX. The muon trigger upgrade, with an online mo- 
mentum measurement, will improve the level-1 trigger rejection for both single 
muons (with ap? cut) and muon pairs (with an invariant mass cut). It will also 
improve the high multiplicity background rejection during the final analysis. 
The muon trigger upgrade is presently under construction. 



2.2 STAR upgrades 



To realize the compelling scientific opportunities in heavy flavor physics, up- 
grades to the STAR detector are required to complete many of the challenging 
measurements. The collaboration has planned a series of upgrades for the near 
and intermediate term to overcome the current shortcomings and enhance its 
heavy flavor capabilities. Implementation of these upgrades will also allow 
optimum utilization of the increased luminosity expected from RHIC II. 



DAQ upgrade and new 
TPC FEE 



Full Barrel Time-of-Flight 
System 



Forward Meson 
Spectrometer (FMS) 




Integrated Tracking * r | 
Upgrade: Active Pixel 
Layers (PXL) & Inner 
Silicon Strip Layers (1ST) 

Fig. 3. Layout of the STAR experiment in 2005/2006 (modified from Ref. [58], 
reprinted with permission from Elsevier). The locations of the planned upgrades 
are indicated by as dashed lines. See text for details. 

The current layout of the STAR detector is depicted in Fig. 3. The medium 
term upgrades to the detector relevant for heavy flavor physics include: a full 
barrel Time-of-Flight detector (ToF) replacing the current ToF patch and the 
Central Trigger Barrel (CTB); new front end electronics for the large Time 
Projection Chamber (TPC); an upgrade to the data acquisition system (DAQ- 
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1000), and a tracking upgrade including a barrel section with two inner layers 
of silicon pixel sensors (PXL) and one layer of silicon strip-pad sensors (1ST) , 
replacing the current Silicon Vertex Tracker (SVT). 

The new ToF system covering the full outer barrel of the TPC is being con- 
structed and installed in STAR over the next two years. The system uses the 
multi-gap resistive plate chamber (MRPC) technology developed at CERN 
and will consist of 3840 MRPC modules with 23,000 channels of readout [59]. 
The modules will cover the TPC outer barrel (|?7| < 1, < < 2ir) and will be 
mounted in 120 trays which will replace the existing CTB scintillation counter 
trays and ToF patch. 

The ToF triples the current momentum range over which tt, K, and p can be 
identified, considerably improving charm meson and baryon reconstruction. 
When the ToF measurement is combined with the TPC dE/dx measurement, 
electrons can be cleanly identified from the lowest momentum measured, ~ 
200 MeV/c, up to a few GeV/c. This capability complements the electromag- 
netic calorimeter which works well for momenta above ~ 2 GeV/c. STAR will 
then be able to reconstruct soft to medium momentum electrons with high 
efficiency and purity to make a comprehensive J/ip measurement. The ToF, 
in conjunction with the electromagnetic calorimeter (EMC), also allows STAR 
to implement a high-level trigger scheme to select J ftp — > e + e~decays in pp 
collisions. 

A series of improvements to the STAR data acquisition system over the past 
several years has brought the recorded event rate capability from the original 
design of 1 Hz to 50 — 100 Hz. To acquire the very large data samples and high 
data rates needed for heavy flavor measurements, a further upgrade has been 
initiated to achieve a minimum 1 kHz recorded event rate which could produce 
data volumes that significantly exceed the capacity for analysis and storage. 
The rare-trigger data sets will especially benefit from the upgrade since the 
pipelined architecture being implemented will virtually eliminate the front end 
dead time, allowing STAR to make full use of rare-event triggers such as that 
designed for the T. 

An increase in readout speed will be achieved by replacing the TPC front end 
electronics (FEE), making use of circuits developed for the ALICE experiment 
at CERN, in conjunction with an upgrade of the STAR DAQ. In addition to 
the increased physics capabilities from the DAQ upgrade, the replacement 
of the TPC FEE, specifically the readout boards (RDO) that collect data 
from the FEE boards, will make space for a future precision tracking chamber 
between the TPC end planes and the endcap calorimeter. Replacing the TPC 
FEE also assures that this system can be maintained for the next decade or 
more. The readouts of the other existing detectors, which will remain in place 
for the RHIC II era, are being adapted to the new high speed DAQ with only 
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minor changes. 




Fig. 4. The proposed geometry for the STAR inner tracking upgrade, the Heavy 
Flavor Tracker. From the inner to the outer radii, the 2 cm beam pipe, the two 
PXL layers, the exoskeleton to strengthen the beam pipe, the one 1ST layers and 
the existing SSD layer are shown [60]. 

In order to address heavy quark energy loss and thermalization, it will be nec- 
essary to cleanly identify open charm. The recent results from both STAR and 
PHENIX on the suppression and flow of non-photonic electrons are intriguing. 
However, without an identified charm sample, the contributions from semilep- 
tonic bottom decays and systematic errors on background subtraction make 
a clear interpretation of these results difficult. Measurement of the yields of 
various charm species will also allow the study of charm hadrochemistry. 

Efficient topological reconstruction of open charm decays requires a tracking 
"point-back" resolution to the primary collision vertex of ~ 50 /im or less. 
Further, the beam pipe and innermost detector layers must be very thin to 
measure the low pt particles, which comprise the bulk of the cross section, 
and thus minimize the systematic errors in extrapolating the measured yield 
to the total yield. A thin beam pipe and inner detector layers are also key 
elements in efficiently vetoing photon-conversion electrons which, along with 
electron identification from the TPC, ToF and electromagnetic calorimeter, 
will measure the soft lepton and dilepton spectra. STAR is thus developing 
a tracking upgrade for the central rapidity region, the Heavy Flavor Tracker 
(HFT). The essential elements under consideration for the HFT are a new thin, 
small-radius beam pipe (0.5 mm thick, 20 mm radius), two layers of thinned 
(50 /im) CMOS pixel detectors at average radii of 2.5 and 8 cm (PXL) and 
one layer of low mass silicon strip-pad sensors (1ST) at 17 cm, see Fig. 4. The 
existing layer of double-sided silicon strip sensors at a radius of 23 cm (SSD) 
will be kept but upgraded to be compliant with the overall higher readout 
speed (DAQ-1000). The two new 1ST layers will connect tracks from the TPC 
and SSD to hits in the pixel layers. These layers will replace the existing 
three layers of silicon drift detector (SVT). It will be necessary to replace the 
SVT since, when RHIC II becomes operational, the SVT will be over 10 years 
old with a readout too slow to be compatible with the upgraded DAQ. Its 
infrastructure (cables and cooling) also adds undesirable mass in the region 



17 



1 < \rj\ < 2. 



3 Projected RHIC yields 



In this section we present some estimates of the quality of the heavy flavor 
measurements that can be achieved at RHIC with increased luminosities and 
the upgraded detectors, reflected in the heavy flavor yields. 









Up to 2008 


2009 


2011 


2013 


Species 


Energy 


Units 


Obtained 


Projected 


Projected 


Projected 


Au+Au 


200 


fib' 1 


380 


610 


1450 


1820 


PP 


200 


pb- 1 


7.0 


14.6 


31.1 


40 


PP 


500 


pb" 1 




36.5 


78 


100 



Table 1 



The anticipated weekly luminosity delivered by RHIC. When estimating rates, this 
delivered luminosity has to be reduced by a factor that accounts for detector up time 
and collision vertex cuts imposed by the detectors. The RHIC projected luminosities 
are maximum values taken from the RHIC Collider Accelerator Division projections 
for 2009-2013 as of February 23, 2008. They reflect the progressive implementation of 
stochastic cooling and other accelerator upgrades and improvements. The numbers 
in the "obtained" column are the best weekly luminosities observed through Run 8. 

Table 1 is a summary of the weekly-integrated, delivered luminosity estimates 
for RHIC spanning the period when the luminosity upgrade is implemented. 
The year 2009 should be taken to mean RHIC Run 9. The weekly-luminosity 
expectations are based on RHIC Collider Accelerator Division guidance as of 
February 23, 2008. The root-mean square (RMS) of the collision diamond is 
assumed to be 20 cm beginning in 2009. We assume that 80% of the RHIC 
beam is in the central bucket (and thus usable by experiments) prior to 2012, 
when the addition of a 56 MHz RF cavity will bring this up to nearly 100%. 

The PHENIX SVTX barrel is expected to be in place starting in 2010. The 
STAR HFT is expected to be in place starting in 2012. In both cases the 
usable collision vertex range will be limited to ±10 cm by the acceptance of 
the silicon upgrade detectors. Because the RHIC diamond RMS length is ±20 
cm this will result in a reduction in effective luminosity for the experiments. 
This is taken into account when calculating yields. 

The quarkonium cross sections are taken from Ref. [61] with an assumed ip' 
to J/ip ratio of 0.14. The charmonia cross sections are reduced by a factor 
of 0.43 in Au±Au interactions, approximately accounting for the suppression 
measured by PHENIX. No T suppression is assumed. 
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3. 1 Projected PHENIX yields 



Table 2 summarizes the projected PHENIX yields per year for critical heavy 
flavor signals for the period 2009-2013 for 12 week physics runs, and also 
includes the yields observed in RHIC runs to date. The estimated yields are 
based on a number of criteria. The detector acceptances are from PHENIX 
simulations. The minimum bias trigger efficiency for hard processes is assumed 
to be 0.75 for pp and 0.92 for Au+Au interactions. Where appropriate, an 
additional, realistic, level-1 trigger efficiency of 0.8 is used. 

Realistic lepton pair reconstruction efficiencies of 0.8 in pp and 0.4 in Au+Au 
collisions are used. An additional efficiency factor of 0.4 is assumed for a 1 
mm displaced vertex cut to identify B — > J/t/jX decays. The PHENIX vertex 
detector is assumed to be in place after 2010, requiring a collision vertex cut 
of ±10 cm. 

3.2 Projected STAR yields 

Table 3 shows a summary of the STAR projected yields per year for vari- 
ous critical heavy flavor signals for RHIC and RHIC II. The detector accep- 
tances and efficiencies are from STAR simulations, or, where available, derived 
from existing measurements. The assumptions of luminosity, interaction ver- 
tices, cross sections and suppression effects are identical to those used for the 
PHENIX projections discussed above. 

The upgrade of the inner tracking system is central to the future STAR heavy 
flavor program. It will provide sufficient tracking resolution to determine the 
heavy flavor meson decay vertices, providing the means to virtually unam- 
biguously identify these probes. However, the complexity and volume of the 
data from this subsystem precludes the use of dedicated triggers that are 
based on the displaced vertex identification. Certain measurements, such as 
D° — > i^ ± 7r =F , can only be made by analyzing large quantities of minimum 
bias data. The substantial increase in readout speed from the DAQ upgrade 
(DAQ1000) can record events with a rate of up to 1 kHz, i.e., 10 times the 
current rate. In the projections discussed here, we assume a more conservative 
rate of 500 Hz. In order to minimize tape costs, STAR plans to filter these 
events offline in computing farms either before or after transfer to the RHIC 
Computing Facility (RCF). The D° yields reflected in Table 3 are therefore 
limited not by luminosity but by bandwidth. 

The same bandwidth limitation holds for the J/ip trigger in A+A collision. The 
STAR J I ip trigger deployed in pp events is based on the coincident signal above 
threshold of two calorimeter towers in conjunction with signals in the referring 
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Table 2 



The projected yields of several heavy flavor signals in PHENIX for 12 week physics 
runs at RHIC. The yields are shown for both pp and Au+Au collisions at y / ijv r /v= 
200 GeV. The approximate yields obtained at RHIC to date are also shown. These 
reflect the fact that RHIC had not yet achieved the full luminosity development for 
Au+Au by Run 4, or for pp by Run 5. The Run 6 pp and Run 7 Au+Au data have 
not yet been fully analyzed. The projected RHIC values for 2011 and 2013 assume 
that the PHENIX SVTX detector is in place, limiting the usable collision vertex 
range to ±10 cm. The collision diamond RMS will be ±20 cm. The reduction in 
usable luminosity is compensated by improvements in signal to background ratio 
for most measurements. The Xc — ► /" + >~7 measurement requires the NCC. It is 
assumed that one NCC arm is installed for 2011 and both are installed for 2013. 

time-of-flight (ToF) segments to veto photon showers. The large multiplicity 
in A + A events, however, renders this scheme impossible. The J/ip A + A 
yields are therefore bandwidth limited while those in pp depend directly on 
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Table 3 



The projected yields of several heavy flavor signals in STAR for a physics run of 12 
weeks at RHIC at y / ijv r /v= 200 GeV. The approximate yields obtained to date are 
also shown. The projected values assume that the STAR Heavy Flavor Tracker, the 
full ToF barrel, and the DAQ upgrades are in place for 2013. From that time on 
the inner tracking limits the usable collision vertex range to ±10 cm. The collision 
diamond RMS will be ±20 cm. Projections for signals marked (min. bias) are based 
on minimum bias events recorded with a rate of 500 Hz and a 33% combined ex- 
periment and machine duty cycle. Dedicated triggers cannot be efficiently deployed 
for these measurements. The corresponding signal-to-background ratios (S/B) are 
also given for the D° measurements. 

the delivered luminosity. However, detailed studies from recent runs indicate 
that a slightly modified J /ip trigger can be effectively deployed for J/^'s with 
p T > 5 — 6 GeV/c. Projections for this specific trigger are shown in Table 3. 

The bottom measurement, B — > J/ipX — > e + e _ , is based on a high-pr J/ift 
trigger, thereby avoiding the issues inherent to the low-p^ trigger. The pro- 
jections shown in Table 3 are based on preliminary results from high-p^ J/tp 
measurements in pp, folding in the expected efficiencies of the inner tracking 
system. The p^-dependent ratio of J/ip from B decays to directly produced 
J /ip was derived from calculations described later. No bottom quark energy 
loss was assumed for the Au±Au yields. 

The decay electrons from T — > e + e~ are at considerably higher pt than those 
of the J/ip, making the calorimeter-based T trigger highly efficient (e > 90%). 
The combination of an efficient trigger and large acceptance makes the study 
of the bottomonium in pp and A + A the main point of the STAR quarkonium 
program and complements PHENIX's strong charmonium program. 
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Table 4 



Projected heavy flavor yields in PHENIX and STAR for a 12 week y/s = 500 GeV 
pp run at RHIC. The projected PHENIX yields for 2011 and 2013 assume that the 
PHENIX SVTX detector is in place, limiting the usable collision vertex range to 
±10 cm. The STAR HFT is assumed to be in place from 2012, and it also will limit 
the usable collision vertex range to ±10 cm. The collision diamond RMS will be 
±20 cm. The Xc M + /U~7 measurement requires the NCC. It is assumed that one 
NCC arm is installed for 2011 and both are installed for 2013. 

3.3 Yields at higher energies 

Table 4 summarizes the expected PHENIX and STAR heavy flavor yields per 
year in pp collisions at y/s = 500 GeV. Although not directly comparable with 
heavy-ion yields from collisions at y / s/v r /v = 200 GeV, the order of magnitude 
larger heavy flavor yields at 500 GeV should help in the understanding of the 
reaction mechanism in pp collisions. 

Table 5 contains a summary of the projected yields per year from the LHC 
detector collaborations for various critical heavy flavor signals in 10 6 s of data 
taking at ^/snn = 5.5 TeV Pb±Pb run, the standard planning number for a 
one year run [62-64]. The ALICE muon spectrometer covers only one side of 
forward rapidity. 

Comparison of Tables 2, 3 and 5 reveal that the projected heavy flavor yields 
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Table 5 



The estimated LHC heavy flavor yields for a 10 s Pb+Pb run at tJsnn = 5.5 TeV 
with 500 /ib _1 integrated luminosity (one year), reported by the LHC experiments. 
As for the RHIC tables, the T rates include all three states. The ALICE yields [62] 
assume binary collision scaling and include shadowing. The ALICE dielectron yield 
estimates are for the 10% most central events while the D° yields are for 10 7 central 
event triggers. The CMS yields [63] assume dN c h/drj = 2500 and include shadowing 
with no quarkonia suppression. For ATLAS [64] the range of the J/ip yield corre- 
sponds to different assumed muon pr trigger thresholds. Good separation of T and 
T' states is expected in ATLAS for |r?| < 2.0. 

for one year of running are similar at the LHC and at RHIC II. The much larger 
heavy flavor cross sections at the higher LHC energy are largely compensated 
at RHIC II by integrated luminosities that result from three times longer runs 
and a factor of 10 higher luminosity. 



4 Open heavy flavor 

In this section we present a more detailed discussion of the theoretical moti- 
vation for studying open heavy flavor in heavy-ion collisions, of the present 
experimental and theoretical status, and of the proposed experimental pro- 
gram of open heavy flavor measurements at RHIC II. 

As described in the introduction, dense matter effects in nuclear collisions may 
change the kinematic distributions and the total cross sections of open heavy 
flavor production. Effects such as energy loss and flow can significantly modify 
the heavy flavor p T distributions but do not, in fact, change the total yields. In 
a finite acceptance detector, however, the measured yields may appear to be 
enhanced or suppressed, depending on the acceptance. Energy loss steepens 
the slope of the heavy flavor p? distribution because the heavy quark p? is 
reduced. If the momentum is reduced sufficiently for the quarks to be stopped 
within the medium, the heavy quarks can take the same velocity as the sur- 
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rounding medium and 'go with the flow'. The present RHIC results on Raa 
and t>2 for heavy flavor decays to leptons show that these effects are indeed 
important for charm quarks and even for bottom quarks. However, higher p? 
measurements and reconstructed charm hadrons are needed to solidify and 
quantify the results. In addition, clean separation of leptons from charm and 
bottom decays is necessary to provide a sensitive measurement of the electron 
Raa from bottom decays. 

Effects that may modify the total heavy flavor yields are the initial parton 
distributions in the nuclei and secondary charm production in the medium. 
The parton distribution functions needed for perturbative QCD calculations 
of heavy flavor production are modified in the nucleus, as was observed in 
nuclear deep-inelastic scattering [65]. At very small momentum fractions, x, 
the gluon fields may be treated as classical color fields. The modifications of 
the parton distributions in nuclei relative to free protons would affect the total 
yields. The effect is expected to be small at midrapidity and moderate pt at 
RHIC but is likely to be more important at large rapidity where lower x values 
are probed. Although thermal charm production from the medium is likely to 
be small at RHIC energies, it could moderately enhance the total yields. 

Since the J/ip yields may be enhanced in nuclear collisions by coalescence 
of uncorrelated c and c quarks in the medium, it is important for charmo- 
nium production in heavy-ion collisions to be properly normalized. The ratio 
of J /if) to open charm production in pp collisions is not a strong function of 
energy. Thus the total charm yield sets the scale against which J/ip suppres- 
sion relative to enhancement can be quantified, as discussed in more detail in 
Section 5. 

4-1 Open Heavy Flavor Theory 
4-1.1 Theoretical Baseline Results 

We now discuss the most recent theoretical baseline calculations of the trans- 
verse momentum distributions of charm and bottom quarks, the charm and 
bottom hadron distributions resulting from fragmentation and, finally, the 
electrons produced in semileptonic decays of the hadrons [6]. 

The theoretical prediction of the electron spectrum includes three main com- 
ponents: the Pt and rapidity distributions of the heavy quark Q in pp colli- 
sions at y/s = 200 GeV, calculated in perturbative QCD; fragmentation of 
the heavy quarks into heavy hadrons, H, described by phenomenological in- 
put extracted from e + e~ data; and the decay of H into electrons according to 
spectra available from other measurements. This cross section is schematically 
written as 
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where the symbol ® denotes a generic convolution. The electron decay spec- 
trum, f(H — > e), includes the branching ratios. 




p T (GeV) p T (GeV) 



Fig. 5. Left-hand side: The theoretical uncertainty bands for c quark and D meson p? 
distributions in pp collisions at = 200 GeV, using BR(c — > D) = 1. The final [66] 
STAR d+Au data (scaled to pp using N co n = 7.5) are also shown. Right-hand 
side: The same for b quarks and B mesons. Modified from Ref. [67], reprinted with 
permission from Elsevier. 

The distribution Ed 3 a(Q)/dpn is evaluated at Fixed-Order plus Next-to- 
Leading-Log (FONLL) level [68]. In addition to including the full fixed-order 
NLO result [69,70], the FONLL calculation also resums [71] large perturba- 
tive terms proportional to a™ \og k (pt / ttiq) to all orders with next-to-leading 
logarithmic (NLL) accuracy (i.e. k = n, n — 1) where m is the heavy quark 
mass. The perturbative parameters are the heavy quark mass, tuq, and the 
value of the strong coupling, a s . Since the FONLL calculation treats the heavy 
quark as an active light flavor at px » m, the number of light flavors used to 
calculate a s includes the heavy quark, i.e. nit + 1 where, for charm, nn = 3 {u, 
d and s). The same number of flavors, n\f + 1, is also used in the fixed-order 
scheme where the quark mass is finite. The central heavy quark masses are 
m c = 1.5 GeV/c 2 and = 4.75 GeV/c 2 . The masses are varied in the range 
1.3 < m c < 1.7 GeV/c 2 for charm and 4.5 < < 5 GeV/c 2 for bottom to esti- 
mate the mass uncertainties. The five-flavor QCD scale is the CTEQ6M value, 
A^ 5 ) = 0.226 GeV. The perturbative calculation also depends on the factor- 
ization (hf) and renormalization scales. The scale sensitivity, a measure 
of the perturbative uncertainty, is calculated using n 2 R F = $ = p\ + m 2 Q as 
the central value while varying \i F and [Mr independently within a 'fiducial' 
region defined by hr^f = £,r,f^o with 0.5 < £n t p < 2 and 0.5 < £,r/£f < 2 
so that {(U£f)} = {(1,1), (2,2), (0.5,0.5), (l,o'.5), (2,1), (0.5,1), (1,2)}. The 
envelope containing the resulting curves defines the uncertainty. The mass and 
scale uncertainties are added in quadrature so that 
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where C is the distribution for the central value, //, max (//, min) is the 
maximum (minimum) cross section obtained by choosing the central value 
with the scale factors in our seven fiducial sets, and m, max (m, min) is the 
maximum (minimum) cross section obtained with £r = £p = 1 and the lower 
and upper limits on the quark mass respectively. 

This range of inputs leads to a FONLL total cc cross section (with n\t +1 active 
flavors) in pp collisions of <t£? nll = 256t^e /ib at ^ = 200 GeV. The the- 
oretical uncertainty is evaluated as described above. The corresponding NLO 
prediction is 2441^ /A>. The predictions in Ref. [15], using m c = 1.2 GeV/c 2 
and hr = fip = 2/io gives cr^ LO = 427 /ib, within the uncertainties. Since 
the FONLL and NLO calculations tend to coincide at small p?, which dom- 
inates the total cross section, the two results are very similar. Thus the two 
calculations are equivalent at the total cross section level, within the large 
perturbative uncertainties. The total cross section for bottom production is 

^FONLL = l g 7l 0.99 ^ 

When the total cc cross section is calculated with the next-to-leading order 
scaling functions of Ref. [69], it does not depend on any kinematic variables, 
only on the quark mass, m, and the renormalization and factorization scales 
with central value [i^p — A*o — m - The heavy quark is always considered 
massive in the calculation of the total cross section and is thus not an active 
flavor in the production calculation. Therefore, the number of light quark 
flavors, nif, does not include the heavy quark. With n\t light flavors and a 
fixed scale, the charm and bottom NLO total cross sections at ^J~s = 200 GeV 
are <7^° nif = 30ll^° fib and ^ LOn » = 2.06^f /ib, respectively [7]. While 
the central values are only about 25% and 10% higher respectively than the 
FONLL results due to a fixed rather than running scale in a s , the uncertainty 
is considerably larger, especially for charm. The larger upper limit on the 
uncertainty is primarily due to the behavior of the CTEQ6M gluon density 
at factorization scales below the minimum scale of the parton density and 
the larger value of a s for n\f light flavors [7]. Other parton densities, with 
lower initial scales and smaller Aq CD give a reduced upper limit. When both 
the NLO and FONLL calculations use the same number of light flavors, the 
results for the total cross section and its uncertainty are in agreement [72]. 

The fragmentation functions, D(c — > D) and D(b — > B), where D and B 
indicate a generic admixture of charm and bottom hadrons, are consistently 
extracted from e + e~ data in the FONLL context [73]. 
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The measured spectra for primary B — > e and D — > e decays are assumed to 
be equal for all bottom and charm hadrons, respectively. The contribution of 
electrons from secondary B decays, B — > D — > e, was obtained by convoluting 
the D — > e spectrum with a parton-model prediction of b — > c decay. The 
resulting electron spectrum is very soft, giving a negligible contribution to the 
total. The decay spectra are normalized using the branching ratios for bottom 
and charm hadron mixtures [74]: BR(B -»• e) = 10.86 ±0.35%, BR(D -»• e) = 
10.3 ± 1.2%, and BR(5 -> D -> e) = 9.6 ± 0.6%. 
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Fig. 6. Left-hand side: The theoretical uncertainty bands for D — ► e (solid), -B — ► e 
(dashed) and -B — > D — > e (dot-dashed) as a function of in y^ = 200 GeV pp 
collisions for \y\ < 0.75. Right-hand side: The final electron uncertainty band in pp 
collisions is compared to the PHENIX [75] and STAR [19, 66] data. Modified from 
Ref. [67], reprinted with permission from Elsevier. 

The left-hand side of Fig. 5 shows the theoretical uncertainty bands for c 
quarks and D mesons, obtained by summing the mass and scale uncertainties 
in quadrature. The band is broader at low px due to the large value of a s 
and the behavior of the CTEQ6M parton densities at low scales as well as 
the increased sensitivity of the cross section to the charm quark mass. The 
rather hard fragmentation function causes the D meson and c quark bands to 
separate only at pt > 9 GeV / c. The right-hand side of Fig. 5 shows the same 
results for b quarks and B mesons. The harder b —>■ B fragmentation function 
causes the two bands to partially overlap until pt — 20 GeV/c. 

Figure 6 shows the individual uncertainty bands for the D — > e, B — > e and 
B — » D — > e decays to electrons on the left-hand side and compares the RHIC 
data to the total band on the right-hand side. The upper and lower limits of 
the total band are obtained by summing the upper and lower limits of each 
component. The secondary B — » D — > e spectrum is extremely soft, only 
exceeding the primary B — > e decays at pt < 1 GeV / c. It is always negligible 
with respect to the total yield. While, for the central parameter sets, the B — ► e 
decays begin to dominate the D — >• e decays at px — 4 GeV/c, a comparison 
of the bands shows that the crossover may occur over a rather broad range of 
electron pt, assuming that the two bands are uncorrelated. 
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4-2 Models of Heavy Quark Energy Loss 

While the heavy quarks are in the medium, they can undergo energy loss by 
two means: elastic collisions with light partons in the system (collisional) and 
gluon bremsstrahlung (radiative). We will briefly review some of the predicted 
results for —dE/dx of heavy quarks for both collisional and radiative loss. We 
then show the predicted effect on the charm and bottom contributions to 
single electrons at RHIC [10]. 

The collisional energy loss of heavy quarks through processes such as Qg — > 
Qg and Qq — > Qq depends logarithmically on the extremes of the heavy 
quark momentum, —dE/dx oc ln(g max /g min ). Treatments of the collisional loss 
vary with the values assumed or calculated for the cutoffs. These cutoffs are 
sensitive to the energy of the heavy quark and the temperature and strong 
coupling constant in the medium. Thus the quoted value of the energy loss is 
usually for a certain energy and temperature. The calculation was first done 
by Bjorken [76] who found —dE/dx 0.2 GeV/fm for a 20 GeV/c quark at 
T = 250 MeV. Further work refined the calculations of the cutoffs [77-79], 
with similar results. Braaten and Thoma calculated the collisional loss in the 
limits E <C mh/T and E ^> itlq/T in the hard thermal loop approximation, 
removing the cutoff ambiguities. They obtained —dE/dx ~ 0.3 GeV/fm for 
a 20 GeV/c charm quark and 0.15 GeV/fm for a 20 GeV bottom quark at 
T = 250 MeV [80]. 

Other models of heavy quark energy loss were presented in the context of 
J/ if) suppression: Could a produced cc pair stay together in the medium long 
enough to form a J/ijrt Svetitsky [81] calculated the effects of diffusion and 
drag on the cc pair in the Boltzmann approach and found a strong effect. 
The drag 2 stopped the cc pair after traveling about 1 fm but Brownian dif- 
fusion drove the quarks apart quickly. The diffusion effect increased at later 
times. Svetitsky essentially predicted that the heavy quarks would be stopped 
and then go with the flow. His later calculations of D meson breakup and 
rehadronization [16] while moving through plasma droplets reached a similar 
conclusion. Koike and Matsui calculated energy loss of a color dipole mov- 
ing through a plasma using kinetic theory and found —dE/dx ~ 0.4 — 1.0 
GeV/fm for a 10 GeV/c QQ [82]. The collisional loss was thus predicted to 
be small, less than 1 GeV/fm for reasonable assumptions of the temperature. 
The loss increases with energy and temperature. Using the hard thermal loop 
approach, Mustafa et al. found —dE/dx « 1 — 2 GeV/fm for a 20 GeV/c 
quark at T = 500 MeV [83]. 

The first application of radiative loss to heavy quarks was perhaps by Mustafa 

2 The drag coefficient A(p 2 ) is related to the energy loss per unit length by A(p 2 ) = 
(-dE/dx)/p 2 . 
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et al. [83]. They included the effects of only a single scattering/gluon emission, 
Qq Qqg or Qg — > Qgg. In this case, the loss grows as the square of the 
momentum logarithm, ln 2 (g max /g min ), one power more than the collisional loss, 
but is of the same order in the strong coupling constant [80] . Thus the radiative 
loss is guaranteed to be larger than the collisional in this approximation. The 
heavy quark mass enters these expressions only in the definition of g max so 
that the mass dependence of the energy loss is rather weak. They found, for 
a 20 GeV/c quark at T = 500 MeV, -dE/dx « 12 GeV/fm for charm and 10 
GeV/fm for bottom. 

These large values suggested that energy loss could be quite important for 
heavy quarks. If true, there would be a strong effect on the QQ contribution 
to the dilepton continuum. Shuryak [84] was the first to consider this pos- 
sibility for A + A collisions. He assumed that low mass QQ pairs would be 
stopped in the medium, substantially suppressing the dilepton contribution 
from these decays. Heavy quarks are piled up at low and midrapidity if 
stopped completely. However, stopped heavy quarks should at least expand 
with the medium rather than coming to rest, as discussed by Svetitsky [81]. 
Lin et al. then calculated the effects of energy loss at RHIC, including ther- 
mal fluctuations, for constant —dE/dx = 0.5 — 2 GeV/fm [14]. These results 
showed that heavy quark contributions to the dilepton continuum would be 
reduced but not completely suppressed. In any case, the energy loss does not 
affect the total cross section. 

Dokshitzer and Kharzeev pointed out that soft gluon radiation from heavy 
quarks is suppressed at angles less than # = m Q/E [13]. Thus bremsstrahlung 
is suppressed for heavy quarks relative to light quarks by the factor (1 + 
0q/6 2 )~ 2 , the 'dead cone' phenomenon. The radiative energy loss of heavy 
quarks could then be quite small. However, Armesto et al. [12] later showed 
that medium-induced gluon radiation could 'fill the dead cone', leading to 
non-negligible energy loss for heavy flavors. They also found that the energy 
loss would be larger for charm than bottom quarks. 

So far the RHIC heavy-ion measurements are not for heavy flavored hadrons 
but for the electrons from their semileptonic decays. If the effects of energy loss 
are substantially different for charm and bottom quarks, then the results in 
Fig. 6 which show that, at high p?, the single electron spectrum is dominated 
by b decays, would suggest that, if charm quarks lose more energy than bottom 
quarks, 6-quark dominance of the electron spectra would begin at smaller 
values of electron p T in A + A collisions. This would, in turn, limit the electron 
suppression factor, Raa, at moderate pt since the large bottom contribution 
would make Raa larger than expected if the spectrum arose primarily from 
charm quark decays. 

The left-hand side of Fig. 7 is an example calculation comparing the c and 
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Fig. 7. Left-hand side: The differential cross section (per nucleon pair) of charm 
and bottom quarks calculated to NLO in QCD [6] compared to single electron 
distributions calculated with the fragmentation and decay scheme of Ref. [6]. The 
solid, dotted and long dashed curves show the effects of heavy quark energy loss 
with initial gluon rapidity densities of dN g /dy = 0,1000, and 3500, respectively. 
Right-hand side: The ratio of charm to bottom decays to electrons obtained by 
varying the quark masses and scale factors. The effect of changing the Peterson 
function parameters from e c = 0.06, e& = 0.006 (lower band) to e c = = 10~ 5 
(upper band) is also illustrated for correlated b and c scales. From Ref. [10], reprinted 
with permission from Elsevier. 

b distributions at midrapidity, as well as their contributions to single elec- 
trons. Electrons from bottom decays dominate the single electron spectra at 
Pt ~ 5 GeV/c for all gluon rapidity densities. This conclusion is further sup- 
ported by the right-hand side of Fig. 7 where the ratio of charm relative to 
bottom decays to electrons is shown. The crossover region here is rather nar- 
row because fip and hr for c and b are are assumed to be correlated while 
they are uncorrelated in Fig. 6. In all cases, the bottom contribution to single 
electrons is large and cannot be neglected in the computation of single elec- 
tron suppression, shown in Fig. 8. Since bottom energy loss is greatly reduced 
relative to charm [10], the possible effect on the electron spectrum is reduced, 
leading to Raa(pt < 6 GeV/c; e) > 0.5 ±0.1. A calculation by Armesto et al. 
[37], with a somewhat different model of energy loss, showed similar results to 
those in Fig. 8. 

Djordjevic and Heinz found that finite-size effects on radiative energy loss in 
a dynamical QCD medium induce a nonlinear path length dependence of the 
energy loss [85]. The light quark and charm quark loss is similar at RHIC while 
nonlinearities only become large for bottom quarks in LHC conditions. Recent 
calculations by Wicks et al. [86] have revisited the importance of elastic energy 
loss and have shown that this component may make a larger contribution to 
the suppression factor than previously expected. 
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Fig. 8. Single electron attenuation pattern for dN g /dy = 1000, left, and 
dNg/dy = 3500, right. The solid curves employ the fragmentation scheme and lep- 
ton decay parameterizations of Ref. [6] while the dashed curves use the Peterson 
function with e c = 0.06 and e;, = 0.006 and the decay to leptons employed by the 
PYTHIA Monte Carlo. Even for the extreme case on the right, the less quenched 
b quarks dilute Raa so much that the modification of the combined electron yield 
from both c and b decays does not fall below ~ 0.5 — 0.6 near px ~ 5 GeV/c. From 
Ref. [10], reprinted with permission from Elsevier. 

A calculation by Adil and Vitev [87] of collisional dissociation of heavy fla- 
vor hadrons suggests that the D and B meson attenuation in a quark-gluon 
plasma are similar. The B dissociation rate is much faster than the D dis- 
sociation rate due to its shorter formation time so that it can undergo more 
fragmentation/dissociation cycles in the plasma than the D, resulting in a sim- 
ilar Raa for D and B mesons and relatively good agreement with the measured 
nonphotonic electron Raa- The heavy quark flow has not been calculated in 
this approach. 

Moore and Teaney [88] and Rapp et al. [89,90] calculated Raa and the non- 
photonic electron elliptic flow, t>2, in a Langevin model of the time evolution 
of heavy quarks in the medium. Both these groups emphasize that elastic 
(collisional) energy loss should be important at low px relative to radiative 
loss since the boost for heavy flavor hadrons in the medium should not be 
large. Both also find a strong correlation between Raa and t>2. Although the 
approaches differ somewhat, the trends are similar in the two calculations. 

Moore and Teaney [88] calculate the diffusion and drag coefficients for charm 
quarks in perturbative QCD. The diffusion coefficient is proportional to the 
inverse square of the strong coupling constant, a s , e.g. D(2nT) oc a~ 2 . They 
present the effects of a range of values for D(2nT) on Raa and t>2, finding 
the largest effects at high p? for small D(2ttT), corresponding to large a s or 
strong coupling in the plasma. 

Rapp et al. [89, 90] calculated the diffusion and drag coefficients assuming that 
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resonant D and B states in the QGP elastically scatter in the medium. Reso- 
nance scattering reduces the thermalization times for heavy flavors relative to 
those calculated with perturbative QCD matrix elements for fixed a s = 0.4. 
The effect is larger for charm than for the more massive bottom quarks. In- 
cluding these states thus reduces the electron Raa at high p T relative to the 
results in Ref. [10] while increasing the electron v 2 to ~ 10% at p T ~ 2 GeV/c, 
in relative agreement with the data. 

Thus, given sufficiently strong coupling and/or resonant states, both R A a and 
v 2 can be described within transport approaches using elastic scattering. More 
and better data is necessary to distinguish the two approaches. 



4-3 RHIC open heavy flavor measurements to date 

Open heavy flavor production cross sections can be measured by reconstruct- 
ing the invariant mass of the heavy quark hadron from its hadronic decay 
products or by detecting leptons from semileptonic decays of those hadrons. 
While both PHENIX and STAR can measure heavy flavor cross sections by 
either technique, PHENIX has some advantages for semileptonic decay mea- 
surements and STAR has advantages for the hadronic decay measurements. 

In both cases the signal to background ratio can be greatly improved if a 
precise measurement of the decay vertex position is available, since hadrons 
containing c or b quarks typically travel several hundred microns from the 
collision point before decaying. Both PHENIX and STAR have plans to add 
secondary vertex detectors capable of the necessary precision. The PHENIX 
SVTX detectors will be installed starting in 2010, the STAR HFT in 2012. 
In addition to reducing the background rates for open heavy flavor decays 
to leptons and hadrons, the secondary vertex detectors make a clean bottom 
cross section measurement possible using displaced vertex decays to J/ip, given 
sufficient luminosity. 

Open heavy flavor cross section measurements based on semileptonic decays 
of charm and bottom mesons are feasible because a small lepton signal can be 
identified in a very large hadron background. The background lepton sources 
are both small and well enough understood that they can be subtracted to get 
the open heavy flavor signal. However, the loss of information about the decay- 
ing heavy meson due to the recoil kinematics is a disadvantage of semileptonic 
decay measurements. Thus charm and bottom decays cannot easily be distin- 
guished. Open charm measurements using hadronic decay products have two 
advantages: the D meson kinematic properties are reconstructed and sepa- 
ration of charm from bottom is far easier because only a small fraction of 
D mesons arise from bottom decays [6]. A disadvantage of hadronic decay 
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measurements is the huge combinatorial background present in heavy-ion col- 
lisions. 

PHENIX has measured open heavy flavor yields via semileptonic decays to 
electrons at midrapidity (\r)\ < 0.35) using the Ring Imaging Cerenkov de- 
tector and electromagnetic calorimeter for electron identification. At forward 
and backward rapidity (1.2 < \r)\ < 2.2) the two muon spectrometers are used. 
PHENIX results are available for pp at midrapidity [75] and forward rapidity 
[91] as well as for d+Au [92] and Au+Au at midrapidity [20]. No open charm 
cross sections from hadronic decays have yet been reported by PHENIX since 
the small central arm acceptance is a disadvantage for such measurements. 

STAR has measured open heavy flavor yields at midrapidity < 1.0) via 
semileptonic decays using either a combination of the time projection chamber 
(TPC) and time of flight (ToF) for electron identification or a combination of 
the TPC and the electromagnetic calorimeter. The backgrounds that must be 
subtracted are much larger than they are for PHENIX because of the larger 
photon conversion rates in STAR and the lack of a hadron blind electron 
identifier. This is compensated somewhat by the larger acceptance. STAR 
electron results are available for pp, d+Au and Au+Au collisions [19]. STAR 
has also measured open charm yields in the range \n\ < 1.0 through hadronic 
D meson decays for d+Au [66], Au+Au [93], and Cu+Cu [21] collisions. 

Because the charm cross section is much larger than the bottom cross section 
at RHIC and dominates the semileptonic decay spectrum for pt < 2.5 GeV/c, 
the integrated non-photonic lepton cross section is usually assumed to be equal 
to the charm cross section. 

4-3.1 Baseline measurements 

Before any conclusions can be drawn about the hot, dense final state from the 
results for heavy-ion collisions, some baseline information is required. Data 
from pp collisions are needed to establish the underlying cross sections and 
kinematic distributions for open heavy flavor and p + A data are needed to 
isolate effects due to gluon saturation and the intrinsic kx distributions in the 
colliding nuclei. 

Both PHENIX and STAR have measured charm production cross sections at 
midrapidity. These measurements have been extrapolated to all rapidities to 
yield total cross sections. These total cross sections are compared to results 
at other energies and to pQCD calculations in Fig. 9. The results for d+Au 
and Au+Au collisions are divided by the number of binary collisions, N co n, 
to directly compare with the pp results. The STAR data are from combined 
fits to hadronic and semileptonic decay data. The PHENIX data are from 
semileptonic decay measurements only. 
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The STAR values are somewhat higher than those for PHENIX. The total 
charm cross sections differ by about 1.5 times the combined standard devia- 
tions of the two measurements, obtained by adding all statistical and system- 
atic uncertainties in quadrature. 

The NLO calculations and their uncertainties are those of the fixed-order 
calculation with three light flavors for charm, as described in Ref. [7]. While the 
larger uncertainty of the fixed-order calculation shown here can accommodate 
the larger STAR cross section, it is not clear which estimate of the uncertainty 
is best. Obviously the large theoretical uncertainty means that there is little 
predictive power in the total charm cross section. The FONLL estimate of the 
uncertainty in the pt distribution is more reliable for finite px [7]. 
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Fig. 9. Comparison of total cross section measurements. The STAR and PHENIX 
results are given as cross section per binary collisions. Vertical lines reflect the statis- 
tical errors, horizontal bars indicate the systematic uncertainties (where available). 
The NLO calculations and the depicted uncertainty bands are described in the text. 

Since charm and bottom quarks are expected to be produced only in the initial 
nucleon-nucleon interactions, their yield should scale as the number of binary 
collisions, iV co n. The left-hand side of Fig. 10 shows the PHENIX measurement 
of the charm invariant yield in Au+Au collisions, scaled by N co n, at midra- 
pidity as a function of Ap ar t [20]. The PHENIX data integrated over pt > 0.3 
GeV/c are consistent with no N paxt dependence, as expected. The effect of the 
final state medium can be seen when the distributions are integrated above 3 
GeV/c. These pr-dependent effects are discussed in the following section. The 
comparison of charm production measured by STAR in pp, d+Au, Au+Au, 
and Cu+Cu indicates that the total charm cross section scales with N co n, 
as expected for hard processes, see the right-hand side of Fig. 10. The d+Au 
and Au+Au cross sections were further constrained by using the \ow-px muons 
and medium-pr electrons from semileptonic charm decays. The measured cross 
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section for each system and centrality was divided by the appropriate value 
of iV co ii to reflect the pp equivalent cross section [21]. The NLO lines are the 
bounds on the total cc cross section from Fig. 9. 
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Fig. 10. Left: PHENIX measurements of the midrapidity invariant charm yields 
scaled by N co \\ as a function of -/Vp art [20]. More than half of the heavy flavor decay 
electrons included in the integral are above 0.3 GeV/c. Right: STAR measurements 
of the charm cross section derived from identified D° mesons in d+Au, Cu+Cu, and 
Au+Au collisions as a function of iV co ii. 




Fig. 11. Left: Compilation of PHENIX [75] and STAR measurements [19,66] of the 
Pt dependence of the semileptonic decay open heavy flavor cross section from 200 
GeV pp collisions, compared with FONLL calculations [75]. Right: The ratio of the 
data to the FONLL calculation. The band depicts the theoretical uncertainty of the 
calculation. Note, that they are smaller than those from the calculation of the total 
cross-section shown in Fig. 9 and 10 (right) as detailed earlier in this section. From 
Ref. [94]. 



The left-hand side of Fig. 11 shows the measured px dependence of semilep- 
tonic open heavy flavor decays to electrons from 200 GeV pp collisions by 
PHENIX [75] and STAR [19,66]. On the right-hand side, the ratio of the 
measured cross sections to the FONLL calculation [6] is shown. The FONLL 
calculation reproduces the pt slope of the data in both cases. The magnitude of 
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the FONLL cross section is 1.5 times smaller than that of the PHENIX data, 
just at the limit of the theoretical uncertainty band. However, the FONLL 
cross section is a factor of ~ 4 smaller than the STAR data, well outside 
the theoretical uncertainties. The discrepancy of a factor of ~ 2 between the 
PHENIX and STAR cross sections in Fig. 11 is not yet explained. 

A limitation of open heavy flavor measurements using semileptonic decays is 
the difficulty of separating charm and bottom decay contributions. Recently, 
the ratio of charm to bottom contributions has been extracted as a function 
of decay electron p? from pp data by STAR and PHENIX. STAR measure- 
ments were made by fitting measured small angle electron-hadron correlation 
functions to a combination of the D and B correlation shapes given by their 
decay kinematics simulated through PYTHIA. Additional pt integrated mea- 
surements were made using like-sign electron-kaon correlations, exploiting the 
difference in the correlation shape between B and D decays. A new recent mea- 
surement is based on the small angle correlation of semileptonic decay electrons 
with identified D° mesons yielding consistent results with the former method 
[38]. The PHENIX bottom/charm ratio [39] was inferred by fitting the shape 
of the electron-kaon reconstructed mass spectrum with a combination of the 
charm and bottom decay line shapes, taken from PYTHIA simulations. The 
measured STAR and PHENIX b/(c + b) ratios are compared in Fig. 12. They 
are in good agreement and indicate that bottom decays become dominant at 
Pt ~ 3.5 GeV/c. 



Q 

+ 



CD 
CD 



I • e-D° (PYTHIA fit) 


"^pi/? preliminary 


- + e-D° (MC@NLO fit) 


- • e-h, Run5 (PYTHIA fit) 




- A e-h, Run6 (PYTHIA fit) 




~ — FONLL 










p+p Vs=200 GeV 



p+p \ s=200 GeV 



p T (GeV/c) 



■ PHENIX (RUN5&6) y <0 35 

■ PHENIX (RUN6) |y| <0 35 

■ FONLL y=0 
. FONLL 

■ FONLL minimum band y=0 




PH ENIX 

Preliminary 



Electron P,(GeV/c) 



Fig. 12. Measurements of the pt dependence of b/(c + b) determined from semilep- 
tonic heavy flavor decays in 200 GeV pp collisions by STAR [38] (left) and PHENIX 
[39] (right). 



A separate measurement of the bottom and charm cross sections was also 
obtained by PHENIX from fits to the dielectron invariant mass spectrum, 
after removal of all non-heavy flavor contributions [40]. The results are in 
agreement with the cross sections inferred by PHENIX from combining the 
measured bj (c + b) ratios with semileptonic decay electron p T distributions. 
The bottom cross sections determined by PHENIX are summarized in Table 6. 
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Method 


d(j bil d y (^ b ) 


(Mb) 


eh mass shape 
e + e~ mass spectrum 


1.34±0.38lg;^(sys) 


4.61 ± 1.3ll|g(sys) 
3.9 ± 2.5^ (sys) 



Table 6 



Summary of PHENIX [39, 40] measurements of the open bottom cross section in 
200 GeV pp collisions. 

Unlike the J/ip measurements discussed in the next section, the current d+Au 
and pp open heavy flavor results are not precise enough for any conclusions 
to be drawn about either shadowing or broadening. Obtaining more pre- 
cise open heavy flavor baseline results is an important priority for the RHIC 
program over the next few years. 

4-3.2 Heavy-ion measurements 

Both PHENIX and STAR have released striking results on suppression of 
single electrons from open heavy flavor decays in central Au+Au collisions. 
PHENIX also has results for the electron t> 2 from open heavy flavor decays. 

The nuclear modification factors and elliptic flow parameters for electrons 
from semileptonic decays of open heavy flavor in Au+Au central collisions 
from PHENIX [20] are shown in Fig. 13. Figure 14 [19] shows the STAR R AA 
measured in Au+Au and d+Au collisions. The Au+Au Raa data from the two 
experiments are in reasonable agreement. Both show very strong suppression 
in central collisions at high p?. When comparing the non-photonic electron 
Raa data to theory, recall that while the electron data contain contributions 
from both charm and bottom decays, the bottom contribution is expected to 
dominate for p T above ~ 4 GeV/c [6]. Since energy loss is predicted to be 
weaker for bottom, the similarity of the electron Raa, Raa ~ 0.2 — 0.3, to 
light-quark hadrons [4] is surprising. 

In Fig. 13, the PHENIX Raa and t>2 data are compared with models that 
calculate both quantities simultaneously. Curve I is a perturbative QCD cal- 
culation with radiative energy loss [37]. It describes the Raa data using a large 
transport coefficient, q — 14 GeV 2 /fm, that also works well for light hadron 
suppression. However, in this model i>2 arises only from the path length de- 
pendence of energy loss and clearly underpredicts the V2 data. Band II is a 
Langevin-based heavy quark transport calculation [90], including elastic scat- 
tering mediated by resonance excitation, also compared to the STAR data in 
Fig. 14. The best simultaneous description of the Raa and t> 2 data is achieved 
with a small heavy quark relaxation time. The curve labeled III is also from a 
transport calculation [88] where the diffusion and drag coefficients are calcu- 
lated in perturbative QCD. The diffusion coefficients required by the data are 
small in both cases, implying a ratio of viscosity to entropy small enough to 
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Fig. 13. PHENIX measurement of the non-photonic electron Raa (upper panel) and 
V2 (lower panel) in central Au+Au collisions [20]. The theory curves are discussed 
in the text. 

be at or near the conjectured quantum bound [20], consistent with estimates 
of elliptic flow and fluctuation analyses for light-quark hadrons [95, 96]. 

The STAR Raa data are compared to model calculations in Fig. 14. The theory 
curves are from calculations discussed previously. The curves labeled I-IV are 
calculations of the electron spectra from both D and B decays incorporating 
final-state c and b quark energy loss. Curve I [10] uses only DLGV radiative 
energy loss with dN g /dy = 1000, while Curve II [37] employs BDMPS ra- 
diative energy loss (Curve I in Fig. 13). Both calculations predict much less 
suppression than observed. Curves III are from a DLGV-based calculation in- 
cluding both collisional and radiative energy loss [86]. While the addition of 
collisional energy loss substantially increases the suppression over the calcula- 
tion with DGLV radiative loss alone, shown in Curve I, it still underpredicts 
the suppression seen in the data. Curve IV is a Langevin-based calculation 
[90], the center of Band II in Fig. 13. Resonance effects reduce the Raa for 
charm-decay electrons to 0.2, but have a smaller effect on bottom quark en- 
ergy loss, causing Raa to increase at higher pt- Curve V is the same as Curve 
II but with only D decays included. 

There is considerable interest in the behavior of the non-photonic electron 
t> 2 for p T > 2 GeV/c where, as discussed earlier, the bottom contribution is 
expected to become important. New preliminary PHENIX semileptonic decay 
t>2 results were shown at Quark Matter 2008. These show that the V2 remains 
high up to 4.5 GeV/c. Together with the observations from PHENIX and 
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Fig. 14. Measurement of the nuclear modification factor of non-photonic electrons 
by STAR for d+Au and Au+Au collisions at different centralities [19]. The theory 
curves are discussed in the text. 
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Fig. 15. PHENIX preliminary measurement of the non-photonic electron v% in mini- 
mum bias Au+Au collisions [39], providing good statistical significance to 5 GeV/c. 
The band is the van Hees (II) result shown in Fig. 13. 



STAR that the bottom cross section starts to dominate the semileptonic decay 
spectra at about 3.5 GeV/c in pp collisions, the existing Raa data and the 
new preliminary v 2 data strongly suggest that bottom quarks experience large 
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energy loss and significant flow at RHIC. 

In spite of the recent pp measurements of the ratio of bottom to charm with 
p T in the semileptonic decay spectra, the lack of electron separation between 
charm and bottom decays remains a serious limitation. Sensitive measure- 
ments of bottom quark Raa and t>2, in particular, in heavy-ion reactions will 
require clean separation of bottom and charm contributions. The ways in 
which this will be addressed in the future RHIC program are discussed in the 
next section. 

4-4 Proposed open heavy flavor experimental program at RHIC II 

Here we focus on the new open heavy flavor physics that becomes available 
with the the RHIC detector upgrades and the RHIC II luminosity upgrade. 

The measurement of charm mesons through their hadronic decay channels 
relies on the efficient reconstruction of the secondary decay vertices. With the 
Heavy Flavor Tracker (HFT) upgrade the STAR experiment will be able to 
measure D mesons down to low p?. However, to-date no high-level trigger 
scheme exist that would allow to trigger on displaced vertices in heavy-ion 
collisions. Consequently, the direct reconstruction of D mesons does require 
large minimum bias data sets and the subsequent offline reconstruction of these 
events. These studies do not depend on large beam luminosity. The situation 
is different for measurements of open heavy flavor through their semileptonic 
decays where the high-p-r lepton can be used to tag candidate events. Triggers 
for high-p T leptons are already deployed in PHENIX and STAR and will be 
mandatory at RHIC-II luminosities. 

The left side of Fig. 16 shows the reconstructed D° — > Kir channel from STAR 
simulations. Note that the signal cannot be distinguished from the background 
(random combinations of K and 7r tracks) without making topological cuts on 
the pairs. The plot illustrates the need for high pointing resolution of the two 
innermost Si-detector layers of the HFT; tight displaced vertex cuts will elim- 
inate most of the prompt hadron tracks from the combinatorial background. 
The overall reconstruction efficiency, taking into account the acceptance, single 
track efficiency and D° reconstruction efficiency increases from ~ 10% at low 
Pt to ~ 25% at high p^. The efficiency is small at low px because the minimal 
boost from the parent particle allows the pool of daughter track candidates 
to be contaminated by primary tracks, diluting the signal and increasing the 
background. 

The right-hand side of Fig. 16 shows the significance of the STAR D° signal for 
10 8 central Au+Au events as a function of pr- With the upgraded STAR data 
acquisition and TPC readout electronics, such a sample is expected in two 
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Invariant Mass (GeV/c 2 ) Pt (GeV/c) 

Fig. 16. Left: Illustration of D — ► i^7r reconstruction with the HFT in STAR. The 
open symbols depict the invariant Kit mass of all candidate tracks without any 
topological cuts applied. The closed symbols show the final signal after secondary 
vertex reconstruction. The dashed line shows the background fit outside the D° 
mass range. The signal peak is shown on a linear scale in the inset. Right: Expected 
significance of the D° measurement in 200 GeV Au+Au collisions with 100 M central 
events as a function of pt- 

weeks of running. The signal significance is directly related to the precision 
with which one can determine suppression, Raa, and flow, t> 2 . 

Expected v2~| 



n— r 




Fig. 17. Left: Estimate by of the D° V2 precision measured by STAR through 
D — ► Kir decays. The error bars depict the anticipated errors of the measurement 
in 10 s Au+Au minimum bias events. The lower curve illustrates V2 model predic- 
tions if only the light quarks in the D mesons flow. The upper curve illustrates the 
increase in V2 if the charm quark also flows [17]. Right: Estimates [97] of PHENIX 
V2 precision for open charm and bottom decay electrons using a displaced vertex 
measurement, calculated for a delivered luminosity of 30 nb~ 1 . The uncertainties 
on bottom decays at the lowest px are dominated by systematic errors. 

A precision measurement of charm elliptic flow is one of the main goals for 
any future RHIC program. In order to simulate minimum-bias Au+Au mea- 
surements, it is necessary to make assumptions about the magnitude and the 
Pt dependence of the anisotropy of the flow distributions. Models of D° flow 
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with and without charm quark flow [17], represent the two possible extremes 
and thus bracket the expected v 2 range, shown on the left-hand side of Fig. 17. 
With the known efficiency and the anticipated signal-to-background ratios, the 
statistical errors on a proposed measurement can be estimated. The small error 
bars shown (symbols) can be achieved by STAR 10 8 minimum-bias Au+Au 
events. The most important momentum range for determining the thermody- 
namic behavior of the D° is pt < 3 GeV/c since other dynamical effects such 
as jet correlations will become important at higher p T . 

Another way to measure open heavy flavor is through semileptonic decays. 
The azimuthal angle of the D and its decay lepton are closely correlated. 
With the exception of very low pt, the v 2 of the lepton is a good proxy for 
the v 2 of the heavy flavor meson. The disadvantage of this method, however, 
is the difficulty of separating leptons from charm and bottom decays. With 
a displaced vertex measurement and the RHIC II luminosity, separation of 
the charm and bottom contributions to the semileptonic decay spectra can 
be done statistically using the different decay lengths for charm and bottom 
mesons (see the right-hand side of Fig. 19). By analyzing data samples with 
different decay length cuts, the fraction of the signal due to b quarks can be 
varied. The addition of a displaced vertex measurement greatly reduces the 
background for all open heavy flavor measurements, improving both statistical 
and systematic errors. The expected precision for separate v 2 measurements of 
charm and bottom via electron decays at midrapidity in PHENIX are shown 
on the right-hand side of Fig. 17. 

At forward rapidity in PHENIX, displaced vertex measurements can be used 
to measure charm and bottom separately using muon decays. In that case, 
adding a 1 cm displaced vertex cut reduces the muon yields from light hadron 
decays by about one order of magnitude, reducing a major background source 
for single muon measurements, and adding a very tight cut around zero ver- 
tex displacement also greatly reduces the punch-through hadron background. 
Figure 18 shows simulations of the expected precision for charm and bottom 
Raa measurements at forward rapidity in PHENIX in one year of RHIC run- 
ning. The right-hand side of Fig. 17 shows the estimated v 2 precision [97] from 
semileptonic charm and bottom decays in the PHENIX central arms for a run 
with a 30 nb^ 1 delivered luminosity. 

Thus, separate Raa and v 2 measurements as a function of pt and y should be 
possible for both open charm and bottom over a wide rapidity range at RHIC 
II. 

With a displaced vertex measurement and RHIC II luminosity, the B — > J / ipX 
decay channel can provide a very clean measurement of open bottom produc- 
tion by both PHENIX and STAR (see Tables 2 and 3). The displaced vertex 
distributions for prompt J/ip and for B — > J/ipX decays accepted in the 
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Fig. 18. PHENIX estimates of the precision of the heavy flavor semileptonic decay 
Raa obtained using a displaced vertex measurement, calculated for a delivered 
luminosity of 30 nb^ 1 for (left) open charm and (right) open bottom. 
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Fig. 19. Left: Comparison of prompt J/^X displaced vertex distribution (in cm) 
with that from B — > J/ipX decays [55]. Note that the prompt J ftp distribution is 
scaled down by a factor of 100. Right: Decay length distributions (in mm) from 
open charm and bottom simulations [55]. 



PHENIX muon arms are compared on the left-hand side of Fig. 19. The yields 
in Tables 2 and 3 assume a displaced vertex cut of 1 mm. A good measurement 
of the cross section and Raa as a function of pt and y for open bottom produc- 
tion will be possible using B — > J/ipX. Even at RHIC II luminosity, however, 
the yields are not expected to be large enough to permit a t>2 measurement. 
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5 Hidden heavy flavor: quarkonium 



In this section we present a more detailed discussion of the theoretical motiva- 
tion for studying heavy quarkonia in heavy-ion collisions. We also summarize 
the present experimental and theoretical status and describe the proposed 
RHIC II experimental quarkonia program. 

5.1 Theoretical results 

5.1.1 Cross sections in pp collisions 

We discuss quarkonium production in the color evaporation model (CEM) 
which can be used to calculate the total quarkonium cross sections. The CEM 
[98, 99] has enjoyed considerable phenomenological success. In the CEM, the 
quarkonium production cross section is some fraction F c of all QQ pairs below 
the HH threshold where H is the lowest mass heavy-flavor hadron. Thus the 
CEM cross section is simply the QQ production cross section with a cut on 
the pair mass but without any constraints on the color or spin of the final 
state. The produced QQ pair then neutralizes its color by interaction with 
the collision-induced color field — "color evaporation". The Q and the Q either 
combine with light quarks to produce heavy-flavored hadrons or bind with 
each other to form quarkonium. The additional energy needed to produce 
heavy-flavored hadrons when the partonic center-of-mass energy, y/§, is less 
than 2m//, the heavy hadron threshold, is obtained nonperturbatively from the 
color field in the interaction region. Thus the yield of all quarkonium states 
may be only a small fraction of the total QQ cross section below 2m//. At 
leading order, the production cross section of quarkonium state C in an AB 
collision is 



a c = F C J2 J d s j dxidx 2 ff(x 1 ,n 2 ) ff(x 2 ,/J 2 ) a^s) 5(s - xix 2 s) (4) 

where s is the square of the parton-parton center of mass energy, ij = qq or 
gg and &ij(s) is the ij — > QQ subprocess cross section. The total QQ cross 
section takes s — > s in the upper limit of the integral over s in Eq. (4). 

The fraction F c must be universal so that, once it is fixed by data, the quarko- 
nium production ratios should be constant as a function of \/s, y and px- The 
actual value of Fc depends on the heavy quark mass, tjiq, the scale, fi 2 , the 
parton densities, ff-(x,fj, 2 ) and the order of the calculation. It was shown 
in Ref. [61] that the quarkonium production ratios were indeed constant, as 
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expected by the model. 



Of course the leading order calculation in Eq. (4) is insufficient to describe 
high p T quarkonium production since the QQ pair pr is zero at LO. Therefore, 
the CEM was taken to NLO [61, 100] using the exclusive QQ hadroproduction 
code of Ref. [101]. At NLO in the CEM, the process gg — > gQQ is included, 
providing a good description of the quarkonium px distributions at the Teva- 
tron [100]. In the exclusive NLO calculation [101], both the Q and Q variables 
are integrated to obtain the pair distributions. Thus, although ji oc uiq in 
analytic LO calculations, at NLO, /i 2 oc m\ = ttiq + p\ where pr is that of 
the QQ pair, p\ = 0.5(p| Q +Pt^)- 

We use parton densities and parameters that approximately agree with the 
QQ total cross section data, given in Table 7, to determine F c for J/ip and T 
production. The fit parameters [102,103] for the parton densities [104-106], 
quark masses and scales are given in Table 7 while the QQ cross sections 
calculated with these parameters are compared to pp — > QQ and n~p — > QQ 
data in Fig. 20. 
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Table 7 



The parton densities (PDFs), quark masses, thq, and scale to mass ratios, £r = 
n/rriT, used to obtain the 'best' agreement to the QQ cross sections. The quark 
mass is given in GeV/c 2 . The inclusive J/tp production fraction, Fj/^, and the 
inclusive T production fraction, Fx, obtained from the data are also given [50]. 

We now describe the extraction of Fc in Eq. (4) for the individual quarko- 
nium states. This is done using J /ip cross sections measured in pp and p + A 
interactions up to y/s = 63 GeV. The data are of two types: the forward cross 
section, ct(xf > 0), and the cross section at zero rapidity, da/dy\y =0 . All the 
cross sections are inclusive with feed down from Xc and ip' decays. To obtain 
Fj/t], for inclusive J/ip production, the normalization of Eq. (4) is obtained 
from a fit using the cc parameters in Table 7. The comparison of uj/^, to the 
xf > data for all four fits is shown on the left-hand side of Fig. 21. The ratios 
of the direct production cross sections to the inclusive J ftp cross section can 
be determined from data on inclusive cross section ratios and branching frac- 
tions. These direct ratios, Rc, given in Table 8, are multiplied by the inclusive 
fitted Fj/^, also shown in Table 7, to obtain the direct production fractions, 
F^ ir = Fj/^Rq, for the various quarkonium states. 



45 



ji 

(o 
o 

t 

ft 

ft 



10 



3 _ 



10 



,2 _ 



10' 



10 



_- 




~ — 1 — 1 — \ 


1 — — — 


III — | — | — | — | — — | — | — | — pz 

-3. 1 = 


-I : 


U 

r 






^j^'^ - 


I (b) H 


I I I I 


1 I I I 


1 1 1 1 1 1 1 1 1 1 1 



20 30 40 50 60 70 10 

Vs" [GeV] 



15 20 25 

Vs" [GeV] 



30 



t 

ft 
ft 



200 


i i | i i i i | i i i i 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


100 






50 






20 


: (c) 




10 




^ - 


5 






2 


1 1 1 1 1 § 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



-MM 


| 1 1 1 1 | 1 1 1 




1 5 


! 1 yS. 


(d) —- 


i i i i 


1 1 1 1 1 1 1 1 1 


i 1 i i i i 1 i i i i 



100 
50 

20 
10 
5 

2 
1 

100 
50 

20 
10 
5 

2 
1 



J2 

O 
t 
ft 

I 



t 
ft 

I 



20 30 40 50 60 70 15 20 25 30 35 40 

Vs" [GeV] Vs" [GeV] 

Fig. 20. The cc, (a) and (b), and 66, (c) and (d), total cross section data in pp and 
ir~p interactions compared to NLO calculations. In (a) and (b), we show ipl (solid), 
V>2 (dashed), tp3 (dot-dashed) and ipA (dotted). In (c) and (d), we show Tl (solid), 
T2 (dashed), T3 (dot-dashed) and T4 (dotted) [50]. 
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Table 8 

Direct quarkonium production ratios, Rc = t c 1v /a 1 ^? where C = J/ip and T. From 
Ref. [107]. 



The same procedure, albeit somewhat more complicated due to the larger 
number of bottomonium states below the BB threshold, is followed for bot- 
tomonium. For most data below y/s = 100 GeV, the three bottomonium S 
states were either not separated or their sum was reported. No Xi?-integrated 
cross sections were available so that the CEM T cross section were fitted to 
the effective lepton pair cross section at y = for the three T(nS) states. 
The extracted fit fraction, F^xi combined with ay and compared to the data 
for all parameter sets in Table 7, is shown on the right-hand side of Fig. 21. 
Using the individual branching ratios of the T, T' and T" to lepton pairs and 
the total cross sections reported by CDF [108], it is possible to extract the 
inclusive T fit fraction, Fx, given in Table 7. The direct production ratios ob- 
tained in Ref. [36] have been updated in Ref. [107] using recent CDF \b data. 
The resulting direct to inclusive T ratios, Rc, are also given in Table 8. The 
sub-threshold bb cross section is then multiplied by F^ ir = FyR c to obtain 
the direct bottomonium cross sections. 
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Vs~ (GeV) vs~ (GeV) 

Fig. 21. Forward J/-0 (left) and combined y + Y' + Y" inclusive (right) cross sections 
calculated to NLO in the CEM. On the left-hand side, we show ipl (solid), ip2 
(dashed), ip3 (dot-dashed) and ip4 (dotted). On the right-hand side, we show Yl 
(solid), Y2 (dashed), Y3 (dot-dashed) and Y4 (dotted) [50]. 

The total cross sections for the charmonium and bottomonium states in pp 
collisions at y/s = 200 GeV are shown in Tables 9 and 10 respectively. 
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Table 9 

The charmonium cross sections for 200 GeVpp collisions. The inclusive and direct 
J/ip cross sections are both given. 



The energy dependence, shown for both states in Fig. 21, is well reproduced 
by the NLO CEM. All the fits are equivalent for y/s = 100 GeV but differ by 
up to a factor of two at 2 TeV. The high energy Y data seem to agree best 
with the energy dependence of Yl and Y2 although Yl underestimates the 
Tevatron result by a factor of 1.4. 



47 





_inc /„u\ 


cr-Y (nbj 


<7T' (no J 


<7T" (, n L)J 






1 1 


O.OU 


Q A Q 

o.4o 


O 1 Q 

z.lo 


1. 61 




0.04 


T2 


7.54 


3.92 


2.49 


1.51 


8.15 


6.34 


T3 


5.75 


2.99 


1.90 


1.15 


6.21 


4.83 


T4 


4.31 


2.24 


1.42 


0.86 


4.66 


3.62 



Table 10 



The direct bottomonium cross sections for pp collisions at 200 GeV. The production 
fractions for the total T are multiplied by the appropriate ratios determined from 
data. 

5.1.2 Cold nuclear matter effects on quarkonium production at RHIC 

It is essential that the A dependence be understood in cold nuclear matter to 
set a proper baseline for quarkonium suppression in A+A collisions. The NA50 
collaboration has studied the J /if A dependence and attributed its behavior 
to dissociation by nucleons in the final state, referred to as nuclear absorp- 
tion. However, the parton distributions are modified in the nucleus relative to 
free protons. This modification, referred to here as shadowing, is increasingly 
important at higher energies, as emphasized in Ref. [109]. We now discuss the 
interplay of shadowing and absorption in d+Au and A + A collisions at RHIC. 

Shadowing, the modification of the parton densities in the nucleus with re- 
spect to the free nucleon, is taken into account by replacing fj(x, if 1 ) in Eq. (4) 

by Ff(x, fi 2 , b, z) = pa(P, z)S j (A, x, /i 2 , b, z)ff(x, fx 2 ) and adding integrals over 
the spatial coordinates. Here is the shadowing parameterization. The den- 
sity distribution of the deuteron is also included in these calculations but the 
small effects of shadowing in deuterium are ignored. The PHENIX J / ip d+ Au 
data as a function of rapidity show a dependence consistent with nuclear shad- 
owing plus a small absorption cross section, common to all quarkonium states 
[110], of 1 — 3 mb. The J/f> production cross section is calculated in the CEM 
using Eq. (4) with the same mass and scale as in cc production. The calcu- 
lations of the d+Au/ pp and A + A/pp ratios are done at LO to simplify the 
calculations since the LO and NLO ratios are equivalent [8]. 

To implement nuclear absorption in d+Au collisions, the per nucleon produc- 
tion cross section is weighted by the survival probability, [111], 



S abs (6, z) = exp | - J dz'p A (b, z')a ahs (z' - z) j , (5) 

where z is the longitudinal production point and z' is the point at which 
the state is absorbed. The nucleon absorption cross section, <r abs , typically 
depends on where the state is produced in the medium and how far it travels 
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through nuclear matter. If absorption alone is active, i.e. no shadowing so 
that S J ' = 1, then an effective minimum bias A dependence is obtained after 
integrating S ahs over the spatial coordinates. If S ahs = 1 also, cr dA = 2Aa p ^. 
When S ahs ^ 1, cr dA ~ 2A a <r pN where, if a a b s is a constant, independent of the 
production mechanism for a hard sphere nucleus of constant density inside Ra, 
ie. pa = poO(RA-b), a = 1 — 9(X a b s / (167rrQ) with tq = 1.2 fm. The contribution 
to the full A dependence of a from absorption alone is only constant if cr a b s is 
constant and independent of the production mechanism [111]. The observed 
J/ip yield includes feed down from Xc and ip' decays, 




y y y 



Fig. 22. Left-hand side: The J/ip d+Au/pp ratio with EKS98 at 200 GeV as a func- 
tion of rapidity for (a) constant octet, (b) growing octet, (c) singlet, all calculated in 
the CEM and (d) NRQCD. For (a)-(c), the curves are no absorption (solid), = 1 
(dashed), 3 (dot-dashed) and 5 mb (dotted). For (d), we show no absorption (solid), 
1 mb octet/ 1 mb singlet (dashed), 3 mb octet/3 mb singlet (dot-dashed), and 5 mb 
octet/3 mb singlet (dotted). Right-hand side: The J/ip d+Au/pp ratio at 200 GeV 
for a growing octet with <j a b s = 3 mb is compared for four shadowing parameteri- 
zations. We show the EKS98 (solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 
(dotted) results as a function of rapidity. From Ref. [8], reprinted with permission 
from APS. 

The J/ij) may be produced as a color singlet, a color octet or in a combination 
of the two. In color singlet production, the absorption cross section depends 
on the size of the cc pair as it traverses the nucleus, allowing absorption to 
be effective only while the cross section is growing toward its asymptotic size 
inside the target. On the other hand, if the cc is produced as a color octet, 
hadronization will occur only after the pair has traversed the target except 
at very backward rapidity. We have considered a constant octet cross section 
(hadronization outside the medium at all rapidity), as well as a growing octet 
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Fig. 23. The AA/pp ratio with the EKS98 parameterization as a function of y for 
octet (upper) and singlet (lower) absorption. In (a) and (b) we show the Au+Au 
results at 200 GeV while the Cu+Cu results are shown at 200 GeV (c) and (d) 
as well as at 62 GeV (e) and (f). The curves are cr a bs = (solid), 1 (dashed), 
3 (dot-dashed) and 5 mb (dotted). From Ref. [9], reprinted with permission from 
Acta Physica Hungarica. 

cross section which hadronizes to a color singlet inside the medium at backward 
rapidities. For singlets, Sfj^ dir ^ S^ s ^ Sffi s but, with octets, one assumes 
that Sfj^ dir = S^ s = S^y s . If this assumption is relaxed and the octet absorp- 
tion cross sections depend on the final-state size, then af bs > a^ s > o J J^ dir 
and Sjy^ dir > S^° c s > Sffi s . The feed down contributions then effectively lower 
the value of a a b s needed to describe the data. As can be seen in Fig. 22, the 
difference between the constant and growing octet assumptions is quite small 
at large ^/s^tv with only a small singlet effect at y < —2. Singlet absorption is 
also important only at similarly negative rapidities. At other rapidities singlet 
production is dominated by shadowing alone. Finally, we have also considered 
a combination of octet and singlet absorption in the context of the NRQCD 
approach, see Ref. [Ill] for more details. The combination of nonperturbative 
singlet and octet parameters changes the shape of the shadowing ratio slightly. 
Including the singlet contribution weakens the effective absorption. The results 
are shown integrated over impact parameter. The calculations use the EKS98 
shadowing parameterization [112] since it gives good agreement with the trend 
of the PHENIX data. For results with other shadowing parameterizations, see 
Refs. [8,9]. 

Several values of the asymptotic absorption cross section, <7 a b s = 1,3 and 5 mb, 
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corresponding to a = 0.98, 0.95 and 0.92 respectively using Eqs. (5) and (6), 
are shown in Figs. 22 and 23 for d+Au and A + A collisions respectively. These 
values of a a b s are somewhat smaller than those obtained for the sharp sphere 
approximation. The diffuse surface of a real nucleus and the longer range of the 
density distribution result in a smaller value of a a b s than a spherical nucleus. 
As will be seen later, there is good agreement with the trend of the PHENIX 
data [113] for <r a b s = — 3 mb. Work is in progress to quantify the shadowing 
parameterization and absorption cross section more precisely over a range of 
energies [114]. 

The current RHIC data are not sufficiently precise to distinguish between 
J /ip production and absorption in the CEM relative to that in the NRQCD 
approach. However, a measurement of the Xc A dependence may be able to 
clarify the situation [111]. In the CEM, the J/ip and Xc distributions differ only 
in the value of F c . In the NRQCD approach, the J ftp is produced primarily 
in a color octet state while the Xc is produced as a color singlet state. Thus 
while the production of both states would exhibit the same shadowing effect, 
a difference in the J/ip and Xc d+Au/pp ratios due to octet relative to singlet 
absorption may be measurable. 

We now turn to the centrality dependence of J/ip production in d+Au and 
A + A collisions. In central collisions, inhomogeneous (spatially dependent) 
shadowing is stronger than the homogeneous (averaged over centrality) result. 
The stronger the homogeneous shadowing, the larger the inhomogeneity. In 
peripheral collisions, inhomogeneous effects are weaker than the homogeneous 
results but some shadowing is still present. Shadowing persists in peripheral 
collisions in part because the density in a heavy nucleus is large and approx- 
imately constant except close to the surface, and because the deuteron wave 
function has a long tail. Absorption is also expected to be stronger in central 
collisions. 

To study the centrality dependence of shadowing and absorption, the d+Au/pp 
and A + A/pp ratios are presented as a function of N co n, 



where T4 and Tg are the nuclear thickness functions and the inelastic nucleon- 
nucleon cross section, a^ N , is 42 mb at 200 GeV. In Figs. 24 and 25, the N co n 
dependence is shown for several representative rapidities, y — — 2, and 2 
for RHIC. The inhomogeneous shadowing parameterization is chosen to be 
proportional to the path length of the parton through the nucleus [109]. For 
more results, see Refs. [8,9]. 

The dependence of the RHIC ratios on 7V co n is almost linear, as seen in Figs. 24 
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Fig. 24. The ratio d+Au/pp as a function of 7V coll for the EKS98 (a), FGSo (b), 
FGSh (c) and FGS1 (d) shadowing parameterizations. The calculations with EKS98 
and FGSo use the inhomogeneous path length parameterization while that obtained 
by FGS is used with FGSh and FGS1. Results are given for y = —2 (dot-dashed), 
y = (dashed) and y = 2 (solid) at 200 GeV for a growing octet with a a b s = 3 mb. 
From Ref. [8], reprinted with permission from APS. 



and 25. The weakest N co n dependence occurs in the antishadowing region, il- 
lustrated by the y = — 2 result (dot-dashed curve). The overall dependence 
on TVcoii is stronger than that obtained from shadowing alone, described in 
Ref. [109], where inhomogeneous shadowing effects depend strongly on the 
amount of homogeneous shadowing. Relatively large effects at low 
companied by the strongest impact parameter, b, dependence. In the transition 
region around midrapidity at RHIC, the b dependence of the ratio d+An/pp 
due to shadowing is nearly negligible and almost all of the N co u dependence 
at y ~ can be attributed to absorption. The y = — 2 results for color singlet 
production and absorption, in the antishadowing region, are fairly indepen- 
dent of N co \\. The stronger AA/pp effect at y = in Fig. 25(a) and (b) can 
be described by the convolution of antishadowing at y ~ — 2 in i?dAu and 
shadowing at y ~ 2, resulting in the dip in AA/pp at y = shown in Fig. 23. 



52 



1.00 

ft 

^0.75 
< 

^0.50 
0.25 



1 1 1 Mill 1 1 1 1 II 


1 1 Ml 


1 1 1 1 1 1 Mil 


i i i i in i i i i mi i 














1.00 




















"~ — — __ 






v\ 








0.75 


a 










0.50 


< 
< 


(a) Au+Au 




- (b) Cu+Cu 200 GeV _ 


. (c) Cu+Cu 62 GeV _ 






1 


1 


1 1 


1 i 


0.25 





i io ioo 100a 

coll 



10 

N 



100 1 



coll 



10 

N 



100 



coll 



Fig. 25. The ratio AA/pp as a function of N co \\ for a 3 mb octet absorption cross 
section and the EKS98 parameterization at y = (dashed) and y = 2 (solid) for 
Au+Au at 200 GeV (a) and Cu+Cu at 200 GeV (b) and 62 GeV (c). From Ref. [9], 
reprinted with permission from Acta Physica Hungarica. 



5.1.3 Models of quarkonium production in heavy-ion collisions 



In-medium properties of quarkonium from lattice QCD: 

Properties of heavy quarks have been used to characterize "thermal properties 
of the QCD vacuum" ever since the first lattice calculations at non-zero tem- 
perature [115]. Modifications of the interactions between heavy, static quarks 
in a thermal heat bath are clearly reflected by changes of the free energy 
which, in the zero temperature limit, reduces to the heavy quark potential 
[116]. To use this information to analyze thermal modifications of quarko- 
nia requires an intermediate, phenomenological step: the construction of a 
temperature-dependent effective potential which then can be used in a non- 
relativistic Schrodinger equation [117-119] or a more refined coupled-channel 
analysis [120,121]. Quite generically, the potential model analyses suggest a 
sequential suppression pattern where heavy quark bound states dissociate at 
temperatures at which their bound state radii become comparable to the De- 
bye screening radius, illustrated in Fig. 26. Table 11 shows quarkonium disso- 
ciation temperatures from Ref. [122]. 



State 


J/^1S) 


Xc(lP) 


V/(2S) 


T(1S) 


Xb(lP) 


T(25) 


Xb(2P) 


T(35) 


T d /T c 


2.10 


1.16 


1.12 


> 4.10 


< 1.76 


1.60 


1.19 


1.17 



Table 11 



Quarkonium dissociation temperatures [122], illustrating the effects of binding en- 
ergy on the dissociation temperature. 

More recently, the calculation of thermal hadron correlation functions and 
their spectral analysis [124] eliminated some of the ambiguities inherent in 
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Fig. 26. Mean squared charge radii of some charmonium and bottomonium states 
compared to the Debye screening radius, = l/mo, and a related scale, r me d, 
an estimate of the distance beyond which the force between a static QQ pair is 
strongly modified by temperature effects [123]. Open (closed) symbols correspond 
to SU(3) (2-flavor QCD) calculations. From Ref. [116], reprinted with permission 
from Springer- Verlag. 

the potential model approach. The spectral analysis, at least in principle, pro- 
vides an ab-initio approach to the calculation of in-medium properties of heavy 
quark bound states. Its predictive power is reduced only by the application of 
statistical tools like the Maximum Entropy Method (MEM) which, however, 
can be steadily improved with further development of the available comput- 
ing resources and numerical techniques. Predictions based on potential model 
calculations as well as the spectral analysis have been reviewed in recent stud- 
ies that have been performed to analyze prospects for quarkonium studies at 
the LHC [50, 125]. In the following, the analysis of thermal hadron correlation 
functions and the extracted spectral functions are discussed. 



The finite temperature, Euclidean time correlation functions, 
G H (T,r,T) = (J H (T,r)f H (0,d)) , 



(7) 



of hadronic currents, Jh = q{r, r)r#g(T, r), where Th denotes a suitable prod- 
uct of gamma matrices that projects onto the appropriate quantum numbers 
of hadron H, are directly related to the spectral functions, <th(u),T). These 
spectral functions encompass all information about thermal modifications of 
the hadron spectrum in channel H so that 

d3 P _ . * ^ ^.rCOshHr - 1/270) 



G H (r,r,T) = J du —La H (u;,p,T) & 



(2tt) 3 " v ' sinh(cu/2T) 



are directly related to experimental observables. In particular, the spectral 
function in the vector channel, ay (a;, p, T), is directly related to the differential 
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Fig. 27. Spectral functions of t/ c , «//?/> (left) and and XcO, Xd (right) at various 
temperatures below and above T c . The vertical bars give the error on the average 
value of the spectral function in the bin indicated by the horizontal bars. 



cross section for thermal dilepton production, 

dW 5a 2 a v {uj,p,T) 
dud?p ~ 27vr 2 u 2 {e"/ T - 1) ' 



(9) 



Note that the rates obtained using this method do not include any contri- 
butions arising from the feed down of other channels into the vector channel 
[107,126]. 

Some general aspects of the influence of a thermal medium on states with 
different quantum numbers can be deduced from the temperature dependence 
of the thermal correlation functions themselves and does not require the ad- 
ditional step of applying the MEM analysis which is based on probabilistic 
assumptions. Such comparisons show that zero-momentum, thermal hadron 
correlation functions in the ground state channels, i.e. the vector (J /if}, T) and 
pseudoscalar (j] c , channels show little modification in a thermal medium 
up to temperatures T > 1.5 T c . Correlation functions corresponding to radially 
excited charmonium states (xc), however, are already strongly modified close 
to or at T c . 

These generic features are reflected by the spectral functions. Although results 
from different groups currently still differ in detail, there are some general 
trends. The J /if} and r] c remain unaffected by the thermal medium up to 
T = 1.5 T c , shown on the left-hand side of Fig. 27. At higher temperatures it 
is unclear whether the J/tf> already disappears at ~ 1.9 T c [127] or persists as a 
strongly modified resonance up to 2.25 T c [29]. The Xco and Xd both disappear 
at T < 1.1 T c , see the right-hand side of Fig. 27. Finite momentum J/if> states 
show statistically significant but still small modifications for T < 1.5 T c [128] 
due to collisional broadening by higher momentum gluons seen by bound states 
moving relative to the heat bath, see Fig. 28 [129] . Strong J / if) binding above T c 
is also supported by the analysis of spatial correlations [131] and the observed 
insensitivity of the thermal vector and pseudoscalar correlation functions to 
spatial boundary conditions [132]. 
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Fig. 28. The ratio of pseudoscalar correlation functions at non-zero momentum 
above and below T c [128]. Deviations from the horizontal dotted lines for increasing 
p indicate stronger temperature-dependent modifications of the spectral functions. 

Bottomonium studies are considerably more difficult since a larger lattice cut 
off is required to properly resolve these states, particularly for temperatures 
well above 2 T c where the T states are expected to be dissolved. First ex- 
ploratory finite temperature results on bottomonium have been reported for 
temperatures up to ~ 1.5 T c . At this temperature, no thermal modifications 
of the T and r\ h have been observed, as expected. The Xb correlation functions 
are, however, modified at T ~ 1.5 T c , similar to the scalar charmonium case 
at T ~ 1.1 T c . To firmly establish the onset of medium modifications in bot- 
tomonium states, however, requires further refined studies of the bottomonium 
system at lower and higher temperatures. 

More recent calculations of quarkonia screening by Mocsy and Petreckzy [130] 
suggest that the quarkonium states break up at lower temperatures. They 
use 2 + 1 flavor lattice results on color screening to fix the parameters of 
the quarkonium potential used in the calculation of the quarkonium correla- 
tors. The most extreme potential compatible with the lattice results is used 
to derive upper bounds on the quarkonium dissociation temperatures in the 
quark-gluon plasma. In this case, although resonance structures appear in the 
spectral functions, the binding energy drops while the thermal width increases. 
When the width is larger than the binding energy, the state can no longer ef- 
fectively be observed. These calculations suggest that the \c and if?' break up 
at or below T c while the J/ip and higher T family resonances dissolve around 
1.2T C . Only the T(IS') state would persist to ~ 2T C . Further calculations are 
needed to resolve the issue. 



Dynamical Coalescence: 



The production of multiple cc pairs in a single collision introduces a new 
charmonium production mechanism [23]. Charmonium states can be formed 
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in-medium by coalescence of a c and a c quark from independently-produced 
cc pairs. 

In the plasma phase, there are two basic approaches: statistical and dynami- 
cal coalescence. Both these approaches depend on being able to measure the 
quarkonium rate relative to total QQ production. The first calculations in 
the statistical approach assumed an equilibrated fireball in a grand canonical 
ensemble [31,32]. This approach could be reasonable at the high energies of 
the LHC, where the number of produced cc pairs is large. But, at lower ener- 
gies, charm conservation is required since a cc pair is not produced in every 
event. More recent calculations assumed a canonical ensemble only for charm 
production [133-135] . Dynamical coalescence models assume that some of the 
produced QQ pairs which would not otherwise do so can also form quarko- 
nium. This coalescence can take place in the QGP [23, 46] or at hadronization 
[33]. The model includes the rapidity differences, |Ay|, between the Q and Q 
and shows that the enhancement is smaller for large |Ay|. The impact pa- 
rameter dependence of statistical and dynamical coalescence is quite different. 
Statistical coalescence gives the largest enhancement in peripheral collisions 
where the volume of the plasma is small, with only a minor enhancement in 
central collisions. Dynamical coalescence produces a larger enhancement in 
central collisions where the number of QQ pairs per event is greatest but still 
produces a significant effect in peripheral collisions [136]. 

Much smaller enhancements are predicted for secondary quarkonium produc- 
tion in the hadron gas, particularly for the J ftp where the additional pro- 
duction is either small (20 — 60%) [34] or about a factor of two [35] at LHC 
energies and smaller still for RHIC. Larger enhancements may be expected for 
the ip' [34]. The predictions depend strongly on the J/ip + 7r(p) cross sections, 
typically not more than 1 — 2 mb [137]. 

These secondary production models are already testable at RHIC where fac- 
tors of 2 — 3 enhancement are expected from coalescence [23, 135]. Hard scat- 
tering of produced particles from different interactions is related to the idea of 
crosstalk between unrelated interactions [138]. Important crosstalk effects were 
predicted in e + e~ collisions at LEP [138] but were not observed. If secondary 
quarkonium production is found, it would indicate the relevance of such effects. 
Secondary quarkonium may be separated from primary quarkonium, which is 
subject to suppression, by kinematic cuts that take advantage of the fact that 
different mechanisms dominate production in different kinematic regimes. 

Predictions of J/ip production by dynamical coalescence suffer from substan- 
tial uncertainties due to the dependence on the charm quark distributions in 
the medium. In fact, it is possible to turn this uncertainty into an advan- 
tage and probe the medium properties using the observed J/ip momentum 
distributions. Two extremes can be considered [22]: either the initial c quark 
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distributions are unchanged by the presence of the medium or the c quarks 
are thermalized. If the charm quark distributions in the medium are identical 
to those of the initial production process, the interactions of charm quarks 
with the medium would be very weak. In this case, both the J / ip rapidity and 
Pt distributions will be narrower than if no plasma is formed simply because 
the center of mass energy of secondary J/ip production is lower than that of 
the initial nucleon-nucleon interactions. The lower energy results in a reduced 
(p T ) and a narrower rapidity distribution. Thus, instead of the transverse mo- 
mentum broadening expected from initial-state multiple scattering going from 
pp to p + A to A + A, the average p\ in A + A would no longer exhibit the 
monotonic increase seen in pp and p + A interactions for increasing A. On 
the other hand, if the charm quarks are assumed to be in thermal equilibrium 
with the surrounding medium, the charm interaction with the medium would 
be very strong. Any J/^'s produced from thermalized charm quarks flowing 
with the medium would have a Pt distribution with a slope characteristic 
of the temperature of the system at the time they were formed, resulting in 
even narrower rapidity and pt distributions. In either case, the effect would 
be largest in central collisions, reverting to "normal" broadening in peripheral 
collisions where on the order of one or fewer cc pairs will be produced since 
the number of cc pairs scales approximately with the number of collisions. 

The (p T ) of coalescence and suppression has been calculated by several groups 
[48,49, 139]. No realistic calculations of the rapidity distribution from models 
including J/ip coalescence are available for heavy-ion collisions so far but, 
based on the assumption of an underlying open charm distribution peaked at 
y — 0, it is predicted [22] that a strong charm coalescence contribution to J/ip 
production will lead to a narrowing of the rapidity distribution, similar to the 
Pt distribution. 

In order to extract the medium properties from secondary J/ip production, 
systematic studies of J / ip production in pp, p + A and A + A interactions are 
necessary. The pp data determine the intrinsic transverse momentum scale for 
a particular energy while the p + A results fix the level of broadening due to 
cold nuclear matter effects which would then apply to A + A interactions at 
the same energy. 

Models of coalescence, of course, also include J/ip suppression. In addition to 
the screening effects discussed previously, the J/ip can scatter with quarks and 
gluons in the plasma which may break it up more efficiently than screening 
effects alone, especially if temperatures significantly above T c are needed for 
screening to dissociate the directly produced J/ip, as discussed in Ref. [33]. 
At low temperatures, relevant for SPS energies, gJ/ip — > cc is effective for 
J j ip breakup by a thermal gluon. However, at higher temperatures where the 
J/ip should be more loosely bound, inelastic parton scattering, g(q,~q)J/ip — > 
g(q,q)cc, calculated using the leading order matrix elements for gc and gq 
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scattering, is more effective. 
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Fig. 29. Excitation function of the ratio of produced J /if) to the number of cc pairs 
in central heavy-ion collisions for N vaxt = 360. From Ref. [140], reprinted with 
permission from Elsevier. 

The relative importance of J/ip suppression and coalescence will change as a 
function of energy, as shown in Fig. 29 for central collisions of heavy nuclei 
with iVp ar t = 360, from Ref. [140]. The J/ip yield is dominated by primordial 
production at SPS energies and by coalescence at the full RHIC energy. 



Other Models: 



Comovers: The comover model, proposed in Refs. [141-143] and further devel- 
oped by Capella et al, see most recently Ref. [144], describes J/ip suppression 
in nucleus-nucleus collisions by final-state interactions with the dense medium 
created in the collision, the "comovers". The SPS J/ip data are reproduced 
with a small effective dissociation cross section, a co = 0.65 mb [145]. At RHIC 
energies, comover interactions alone would suggest stronger suppression than 
at the CERN SPS since the density of produced particles is larger at RHIC. 
Since the produced particle density is largest at midrapidity, comover inter- 
actions alone would predict stronger suppression at y — than at forward 
rapidity [146]. 

However, although the comover model is based on rate equations, the gain 
term was only previously considered in Ref. [142] and reintroduced by Capella 
et al. in an attempt to account for the stronger observed suppression at for- 
ward rapidity at RHIC [144]. The open charm cross section at forward rapidity, 
where RHIC data are not yet very precise, was obtained from the shape of 
PYTHIA simulations of the rapidity distribution. As expected, coalescence ef- 
fects are stronger at midrapidity so that the behavior calculated in Ref. [146] 
is reversed in Ref. [144] and forward suppression is found to be larger. 
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Sequential Suppression: As discussed previously, interpretation of the quarko- 
nium dissociation temperatures in Table 11 suggests a sequential suppression 
pattern in which the quarkonium states dissolve when the system is at or above 
T d . Applied to RHIC, the ip' and Xc would be completely suppressed above a 
critical energy density, e(T d ) while the temperature at RHIC is too high for 
direct J/ip suppression by color screening. Since the turn on of suppression 
at e(Td) would result in a step-wise suppression pattern, the suppression is 
typically smeared over a range of impact parameters. 

Several models of sequential suppression have been proposed [147-150]. Most 
include some level of normal cold nuclear matter effects but neglect gluon 
dissociation of final-state J/^'s. The (p?) of J/ip suppression, is typically flat 
or decreasing with centrality in these scenarios. 

In Ref. [147], the normal absorption cross section is fixed by the RHIC d+Au 
data in each rapidity region and then applied to the Au+Au data. The thresh- 
old energy density model, originally proposed in Ref. [148], suppresses the J /ip 
yield above a threshold energy density, neglecting Xc and ip' feed down. It was 
recently applied to RHIC but included no rapidity dependence [149] . The most 
dynamically complete model of sequential suppression describes suppression 
in a dynamically-expanding quark-gluon fluid using relativistic hydrodynam- 
ics [150]. 

Conformal Field Theory: Recently, M = 4 super Yang-Mills theory, a conformally- 
invariant field theory, has been applied to J/tp suppression [151]. The screening 
length is described by the dynamics of a color single QQ dipole of velocity v 
moving through the hot plasma. In its rest frame, the dipole feels the medium 
as a hot wind moving past it. The screening length in this medium decreases 
as 7 -1 / 2 or (1 — i> 2 ) -1 / 4 so that the quarkonium dissociation temperature would 
decrease with increasing v {pt)- Thus high p? J/V's could be suppressed by 
the plasma while low pt J/ip's are unaffected. Such novel suppression patterns 
are only observable at RHIC for sufficiently high luminosities [151]. 



5.2 Status of Quarkonium Physics at the CERN SPS 

The prospects of a "clean" QGP signature, destruction of the J/ip by color 
screening, was discussed in the landmark paper by Matsui and Satz in 1986 
[27]. This triggered an extensive experimental program at the CERN SPS. 
HELIOS-III [152] and NA38 [153] (subsequently NA50 [154] and currently 
NA60 [155]) made detailed measurements of the dimuon invariant mass spec- 
trum around midrapidity. Despite early enthusiasm and enormous statistics 
(see Fig. 30), the picture that evolved is still rather ambiguous. The SPS mea- 
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surements must also be understood in light of the many results on quarkonium 
production in p + A collisions from fixed target experiments. The status of the 
SPS program is summarized in this section. 




Fig. 30. Dimuon invariant mass spectrum from 158 ^4GeV Pb+Pb collisions at 
NA50. Modified from Ref. [26], reprinted with permission from Elsevier. 

Feed down contributions (see Fig. 31) from higher charmonium states, Xc — * 
J/ipj (~ 30%) and ip' — > J/ipn-n (~ 10%), are important [156,157]. The Xc 
has not yet been measured by the heavy-ion detectors at the SPS although 
it has been seen in pp and p + A experiments there. These measurements are 
extremely difficult and the large scatter of available data depicted in Fig. 32 in- 
dicates that better measurements are desperately needed. Although the NA60 
experiment is planning to conduct this analysis, the feasibility with the present 
data sets still has to be verified. 



f(2S) 




hadrons hadrons 7* radiative 
jPC = o-+ 1— 0++ 1++ 1+- 2++ 

Fig. 31. Charmonium mass levels and spin states. Common feed down channels are 
indicated. 

The J/ip and ip' have substantially modified cross sections when produced in 
normal (cold) nuclear matter. Considerable effort has gone into studying p + A 
data from six targets covering a wide range of nuclei, from Be to Pb, measured 
with proton beam energies of 400 and 450 GeV [159] . The p + A J/ip and i/j' 
data were fitted assuming a single effective absorption cross section for each 
charmonium state [it. no explicit treatment of feed down effects), and without 
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Fig. 32. The fraction, Ry C i °f observed J/ip's originating from radiative 
Xd,2 J 1^1 decays as a function of energy for proton and pion beams. From 
Ref. [156], reprinted with permission from IOP. 

shadowing (analyses including shadowing effects are now underway [160]) or 
comover absorption [45]. It was concluded that absorption cross sections of 
4.2 ± 0.5 mb and 7.7 ± 0.9 mb for the J/i[) and ip' respectively should be used 
to predict the cold nuclear matter baseline suppression for Pb+Pb collisions 
at 158A GeV beam energy. 

Studies of the A dependence in cold matter were also made at Fermilab over 
a larger xp range, albeit at higher ^/s/vw [161]. A very strong xp dependence 
was observed for xp > 0.2, as depicted in Fig. 33. Effects such as shadowing, 
absorption and energy loss play roles at different xp, resulting in the observed 
dependence [8]. 




Fig. 33. Measurement of J/ip and if)' absorption in 800 GeV p + A collisions as a 
function of xf at the Tevatron [161], reprinted with permission from APS. 

Figure 34 summarizes the SPS results on J/ip and ip' suppression in Pb+Pb 
S+U and p + A collisions [162, 163]. The figure shows measured J/ip and i/j' 
yields (normalized to the Drell-Yan dimuon yield since it is proportional to 
the number of collisions) as a function of the average path length of the cc 
pair through the nucleus, L, a measure of centrality. The measured yields 
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are divided by those expected from cold nuclear matter effects, calculated 
using the effective nuclear absorption cross sections, fit to the p + A data, 
mentioned earlier. The p + A and S+U data in Fig. 34 have been rescaled to 
^s NN = 1 7.3 for comparison with the Pb+Pb data. Isospin corrections were 
applied to the Drell-Yan cross sections. Strikingly, the heavy-ion J/ip data 
show no suppression beyond that expected from cold nuclear matter up to 
L ~ 7 fm, after which suppression beyond cold nuclear matter expectations 
becomes significant. 

Alternatives to QGP models of suppression, such as the comover model, are 
able to describe the observed J/ip suppression by assuming breakup of the 
bound state by comoving matter [45]. Although these approaches make some 
unrealistic assumptions about the hadron density, it is possible that some 
fraction of the observed suppression is due to comover absorption. 

The ip' suppression in Fig. 34 is similar to that of the J / ip although suppression 
beyond cold nuclear matter effects becomes significant at a much shorter path 
length, L ~ 4 fm, e.g. in peripheral S+U collisions. Thus the ip' appears to 
be more readily suppressed in the final state. This may be because the ip', 
lying 50 MeV below the DD threshold, can be more easily broken up by 
interactions in the medium. The strong ip' suppression measured by NA50 has 
been interpreted as both total suppression of the ip' by color screening [164] 
and a larger interaction cross section for comovers [165] . 
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Fig. 34. The ratio of measured J/ip and tp' yields to those expected from cold 
nuclear matter effects in heavy- ion and p + A collisions as a function of the cc path 
length through the colliding nuclei [162]. The measured and expected yields are 
both normalized to the Drell-Yan yields in the mass range 4.2 < M < 7 GeV. The 
lighter systems are scaled to ^/s^ = 17.3 GeV for comparison with the Pb+Pb 
data and isospin corrections are applied to the Drell-Yan yields. 

The NA60 In+In J ftp measurements exhibit a level of suppression similar to 
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the NA50 Pb+Pb data [166]. The J /if) suppression observed in A + A interac- 
tions at the SPS can, for the most part, be accounted for by the assumption 
that the more loosely bound ip' and Xc states are both suppressed by plasma 
production, eliminating their contribution to the inclusive J /if) measurement. 
The direct J/ip contribution is assumed to not be suppressed at the SPS 
[164,167]. However this indirect J /if) suppression picture does not seem to 
be easily reconciled with the strong observed if)' suppression for path lengths 
corresponding to rather peripheral collisions, see Ref. [163], well before any 
reduction in the J /if) yield is observed. 

Charm production was not measured in heavy-ion experiments at the SPS 
until very recently. Figure 30 shows that the open charm contribution to the 
dilepton continuum in the J /if) mass region is negligible at the SPS. Open 
charm measurements are, however, key to understanding the intermediate 
mass dilepton region. The NA60 experiment has used displaced vertices to 
separate charm decays from prompt dileptons. They have presented prelimi- 
nary In+In results which show that while the enhancement in the intermediate 
mass region is confirmed, it is inconsistent with enhanced open charm. Instead, 
the enhancement is consistent with a prompt dilepton source [168]. 

The picture emerging from SPS studies is still somewhat ambiguous. In spite 
of the good precision of the heavy ion data for both J /if) and if)', and the effort 
that went into determining the cold nuclear matter baseline from p + A data, 
the centrality dependence of the suppression pattern is insufficient to draw 
unique conclusions. On the other hand, the vast experience gained at the SPS 
can and should be taken into account at RHIC. The main lesson learned is 
that a simple J / if) measurement in A + A collisions as a function of centrality 
is insufficient to draw unique conclusions, even when a cold nuclear matter 
baseline is available. Rather, a systematic and detailed study of all related 
aspects, i.e., a systematic study of open charm, J/ip, if)', and Xc production 
in pp, p + A, and A + A collisions is required. Centrality, rapidity, and A de- 
pendence studies are mandatory so that competing models based on different 
mechanisms can be differentiated. 



5.3 Quarkonium measurements to date at RHIC 

All of the published J/ip results from RHIC to date are from PHENIX. Some 
preliminary J /if) results from STAR were presented at Quark Matter 2008. 
PHENIX measures quarkonium yields by reconstructing their invariant mass 
from dilepton decays. Dielectrons are used in the central arms (|?7| < 0.35) and 
dimuons are used in the muon arms (1.2 < \rj\ < 2.2). STAR uses dielectrons 
within the TPC acceptance (\rj\ < 1). 
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PHENIX has reported J/ip results at 200 GeV from pp [44,113,169], d+Au 
[47, 113], Au+Au [42] and Cu+Cu collisions [43]. STAR reported preliminary 
high pt J/ip measurements in pp and Cu+Cu collisions at Quark Matter 2008 
[170]. 

PHENIX reported observation of T -> y+ir for I- 2 < \v\ < 2 - 2 in 200 GeV 
pp collisions from RHIC Run 5 [171]. This very low statistics measurement 
(27 counts in both muon arms combined) was used to make a preliminary 
cross section estimate of BR(da/dy) y= i, 7 = 45.2 ± 9.5 (stat) ± 6.3 (sys) nb. 
STAR observed about 50 T — > e + e~ at \y\ < 1, leading to a preliminary 
cross section of BR(da/dy) y=0 = 91 ± 28 (stat) ± 22 (sys) nb [172]. In both 
cases, the cross section is for the lowest three T states combined. It would be 
difficult to make definitive T measurements at present RHIC luminosities, but 
crude measurements of the T nuclear modification factor by both PHENIX 
and STAR are expected to be possible at RHIC with about 10 times the Run 
4 integrated luminosity, which will most likely be available after the next long 
Au+Au run. 



5.3.1 Baseline quarkonium measurements at RHIC 

PHENIX has measured J/ip cross sections in pp [44] and d+Au collisions 
[47, 113] at 200 GeV. The rapidity dependence is summarized in Fig. 35. The 
left side shows the invariant J / ip cross section in pp collisions while the right 
side shows the nuclear modification factor, -RdAu, for minimum bias d+Au 
collisions. 
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Fig. 35. The rapidity dependence of the pp J/ip cross section at 200 GeV [44] (left- 
hand side, reprinted with permission from APS) and the nuclear modification factor 
for minimum bias d+Au collisions [47] (right-hand side). The curves on the left- 
hand side are fits used to extract the total cross section and estimate the systematic 
error. The deuteron is defined to be moving toward positive rapidity. 



The data shown on the right-hand side of Fig. 35 are compared [47] with 
d+Au calculations [8] that include shadowing (EKS98 [112] and nDSg [173]) 
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Fig. 36. The rapidity dependence of the d+Au J /if) invariant yield at 200 GeV 
[47] compared with calculations [8] including shadowing and a J/ if) N breakup cross 
section. The left-hand side compares the data with calculations using the EKS98 
parameterization [112] while the right-hand side employs the nDSg parameterization 
[173]. 



and J/ipN breakup in Fig. 36. The shaded bands show the one standard 
deviation limits of fitting the J/ipN breakup cross section. The extracted 
values are consistent, within large uncertainties, with the 4.2 ± 0.5 mb cross 
section determined at the CERN SPS [159, 174]. 

The J /ip nuclear modification factor in d+Au collisions is shown as a function 
of centrality in Fig. 37 [47] for the forward, central and backward rapidity 
regions covered by the three PHENIX arms. The shaded bands show the same 
calculation as in Fig. 36 but the breakup cross section is now fit to the data 
in each rapidity region. The N co n dependence is the same because the impact 
parameter dependence of the shadowing is assumed to be proportional to the 
J/ip path length through the medium in both cases [8]. It is evident that the 
parameters are not well constrained by the current d+Au data. 



5.3.2 Quarkonium measurements in heavy-ion collisions at RHIC 

Measurements with the statistical precision needed to provide a strong test 
of models of J/ip production in heavy- ion collisions are now available from 
RHIC. These PHENIX measurements are from the Run 4 Au+Au and the 
Run 5 Cu+Cu data sets. The main features are summarized here. 

Figure 38 shows the J/ip nuclear modification factor, Raa, measured in 200 
GeV Au+Au collisions at both central and forward/backward rapidities [42], 
along with the ratio of forward to midrapidity Raa values, Rj^^/R^a - Some 
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Fig. 37. The nuclear modification factor in d+Au collisions measured at forward 
(bottom), central (middle) and backward (top) rapidity as a function of centrality. 
The deuteron is moving toward forward rapidity. The curves, including both shad- 
owing and nuclear breakup, are discussed in the text. From Ref. [47], reprinted with 
permission from APS. 



systematic uncertainties cancel in the ratio. The ratio R^a^/R^aai in the 
lower panel of Fig. 38, shows that when N pauTt > 150 the Au+Au nuclear mod- 
ification factor at forward rapidity is only about 60% of that at midrapidity. 
Thus the Au+Au data show that the suppression is weaker at the maximum 
transverse energy density (midrapidity) . The weaker suppression could be due 
to either entrance channel effects (shadowing) or final-state effects such as J/ ip 
formation by coalescence of uncorrelated charm quarks. 

Figure 39 shows the nuclear modification factor for Au+Au and Cu+Cu at 
central and forward/backward rapidity, along with R f X A waTd / R^a ■ Within un- 
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certainties, the data for Cu+Cu and Au+Au agree where they overlap in iV part . 
It is of considerable interest to see what constraints the d+Au data place on 
the cold nuclear matter contribution to the heavy-ion Raa by using the fits 
discussed in the previous section. However, R t AA aid /^AA * s independent of the 
breakup cross section in those models because the rapidity dependence of the 
absorption model is small in the range of the data [8]. Thus the model calcula- 
tion of i?^™ ard / BF^a reflects only the effects of the assumed shadowing model. 
To explore how the d+Au data constrain the forward to midrapidity ratio for 
heavy ions, an ad hoc model was used to parameterize the d+Au data. In the 
ad hoc model, the breakup cross section is allowed to vary independently at 
y = and \y\ = 1.7. The resulting parameterizations of the d+Au data can 
be used to predict the cold nuclear matter contribution to the heavy-ion Raa 
independently at each rapidity, with independent uncertainties, and allow the 
constraints on R^^/R^f to be quantified. The calculations for EKS98 are 
compared with the heavy-ion data in Fig. 39. The nDSg results are similar. 
The cold nuclear matter Raa predicted by the ad hoc fits to the d+Au data 
are in good agreement with the Cu+Cu Raa measured in more peripheral 
collisions. This is illustrated in Fig. 40, which shows the survival probability 
relative to cold nuclear matter effects on the Cu+Cu data, calculated using the 
ad hoc fits to the d+Au data shown in Fig. 39. Below iV part ~ 50 the survival 
probability is consistent with unity within the ~ 15% uncertainties. Beyond 
that, the uncertainty in the cold nuclear matter reference rapidly worsens. 

The bottom panel of Fig. 39 shows that the existing d+Au data allow the 
additional suppression from hot nuclear matter at forward rapidity to be any- 
where between zero and a factor of two for central Au+Au collisions. Clearly, 
better d+Au data are required before the effects from cold nuclear matter 
are quantified for Au+Au collisions. The Run 8 d+Au data set will provide a 
greatly improved cold nuclear matter reference. 

The left-hand side of Fig. 41 [175] compares the PHENIX Au+Au data to the 
comover model [146], with shadowing and comover interactions, where the 
suppression is predicted to be larger at y = than at y — 1.7, in disagreement 
with the data. The right-hand side of the figures compares the y = Au+Au 
data with calculations containing charmonium dissociation by thermal gluons 
with no coalescence [48, 176]. All of these calculations overestimate midrapid- 
ity J/ip suppression. Figure 42 [175] compares the midrapidity data to model 
calculations [22,48,176-178] including J/ip formation by coalescence. These 
models generally do a better job of describing the magnitude, if not the de- 
tailed centrality dependence, of the suppression. However they rely heavily on 
the rapidity density of charm production as input to the coalescence calcula- 
tions and the charm distributions are currently not very precisely defined by 
the RHIC data, especially at forward rapidity. 

The extended comover interaction model [144], including J ftp coalescence as 
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Fig. 38. The nuclear modification factor as a function of centrality for 200 GeV 
Au+Au collisions measured at central and forward/backward rapidity [42]. 



well as suppression, shown in Fig. 43, has slightly stronger suppression at 
forward rapidity than seen in Fig. 41, in better agreement with the data. A 
midrapidity calculation by Zhao and Rapp [49] extends the suppression and 
coalescence model [140] to include momentum-dependent dissociation rates. 
The J/if)N cross section was reduced to 1.5 mb, closer to the PHENIX d+Au 
data, and a longer charm quark thermalization time, 7 fm/c was assumed to 
better match the magnitude of the observed suppression in the most central 
Au+Au collisions. The calculated centrality dependence of Raa at midrapidity 
is in good agreement with the data. Reference [49] also presents calculations 
of (p T ) an d Raa function of Pt, discussed later. 

PHENIX also showed a measurement of the nuclear modification factor for 
62 GeV Cu+Cu collisions at Quark Matter 2005. While this measurement 
has relatively low statistics, the 62 GeV data exhibit similar, and perhaps 
slightly stronger, suppression in the most central collisions than that seen in 
200 GeV Cu+Cu collisions, consistent with predictions of color screening and 
coalescence [140]. 

There have been recent predictions of cc coalescence effects on the rapidity 
and px dependence of the J /if) yield [48, 49, 139]. Figure 44 shows the average 
J/ip p T , (p T ), obtained from data below 5 GeV/c, as a function of N part for 
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Fig. 39. The nuclear modification factor as a function of centrality for 200 GeV 
Au+Au and Cu+Cu collisions measured at central and forward/backward rapidity 
[42,43]. 



Au+Au and Cu+Cu collisions in the two rapidity regions covered by PHENIX 
[42, 43]. The Au+Au and Cu+Cu data are in very good agreement at the same 
Afpart and rapidity. Also shown are calculations by Yan, Zhuang and Xu [48] 
and Zhao and Rapp [49] that incorporate the effect of coalescence on the 
centrality dependence of {pj). An earlier calculation by Thews [139] predicts 
a larger and more steeply rising J /if) (p\) than shown in Fig. 44. However all 
calculations with J /if) coalescence agree fairly well with the (p\) data, even 
though coalescence is a smaller fraction of the J /if) yield in Refs. [48,49]. 

There has been considerable interest in what might be learned about the J / if) 
production mechanism from Raa{pt), since coalescence formation of the J /if) 
is generally expected to decrease at high p T . PHENIX has published data 
on the pt dependence of Raa in Au+Au and Cu+Cu collisions [42,43]. The 
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Fig. 40. The ratio of measured Raa to the cold nuclear matter Raa calculated from 
the ad hoc fits to d+Au data, shown in Fig. 39, for Cu+Cu collisions. The one 
standard deviation uncertainties in the cold nuclear matter reference are indicated 
by bands around unity. 
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Fig. 41. Left: Nuclear modification factor compared with a comover model cal- 
culation at central and forward/backward rapidity. Note that the model predicts 
maximum suppression at y = 0. Right: Comparison with two calculations of gluon 
dissociation at midrapidity [175]. 



measured Raa for central collisions is flat within errors over the range of the 
data, extending only to 5 GeV/c. Zhao and Rapp [49] have calculated the 
Pt dependence of the J/ip Raa at midrapidity and compared it with the pt 
dependence in Au+Au collisions measured at four different centralities. The 
calculation agrees well with the data, as seen in Fig. 45. Coalescence is most 
important for more central collisions at low pt- It is a negligible effect above 
4 GeV/c. An important goal of future measurements is to improve the poor 
statistical precision of the current data beyond about 3 GeV/c. Any AdS/CFT 
effects [151] would become important for p T > 5 GeV/c, beyond the current 
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Fig. 42. Nuclear modification factor compared with models that include J/tp for- 
mation by coalescence. The calculations are for y = only [175]. See the text for 
details. 




Fig. 43. Calculation of J/i/j Raa at central and forward rapidities in the comover 
model with coalescence [144]. 



range of precision. Extending the Raa measurement with good precision to 
much higher will require RHIC II luminosity, as discussed later. 

A strong charm coalescence contribution to J j ip production will result in a nar- 
rower rapidity distribution, similar to the pt distribution. PHENIX extracted 
the RMS of the rapidity distribution versus centrality from their Au+Au, 
Cu+Cu and pp data [42-44]. The extracted Cu+Cu data are independent of 
centrality within 2 — 3% uncertainties. The extracted Au+Au values decrease 
by 10% with increasing centrality, roughly a two standard deviation reduction. 
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Fig. 44. The (p T ) for 200 GeV Au+Au and Cu+Cu collisions measured by PHENIX 
at midrapidity (left) and forward/backward rapidity (right) [43] compared with 
calculations without coalescence (top) and including coalescence (bottom) [48,49]. 
The pt integrals are restricted to less than 5 GeV/c to reduce the systematic errors. 



In summary, while there is evidence favoring a coalescence contribution to 
the J I ij) yield in central collisions, it does not seem to be definitive. The sup- 
pression as a function of centrality shown in Figs. 42 and 43 is reasonably 
consistent with models that include coalescence, as is the (p T ) in Fig. 44 and 
Raa{pt) shown in Fig. 45. The observed narrowing of the rapidity distribution 
is weak, and its implication will not be clear until the open charm rapidity dis- 
tribution is better determined experimentally and the effects of shadowing and 
other initial state phenomena are understood. The initial-state effects should 
be better constrained by the more precise d+Au data obtained in Run 8. Pre- 
cise open charm rapidity distributions will require upgrades to the PHENIX 
and STAR detectors to identify open charm and bottom decays by displaced 
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Fig. 45. The J/ip Raa at midrapidity as a function of px in Au+Au collisions in 
four centrality bins, compared with calculations including formation by coalescence 
[179]. 



decay vertices. On the other hand, the predicted narrowing of the pt distri- 
bution relative to broadening of the primordial J/xjj px distribution expected 
in heavy-ion collisions is based on a steeply falling charm p? distribution, well 
established by the existing open charm data, and the calculations seem to 
reproduce the data reasonably well. 
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PHENIX is attempting to extract the J /if) t>2 from Run 7 Au+Au data. The 
current restricted heavy-ion data sets place major limitations on the precision 
of a t>2 measurement. A statistically meaningful J/if> polarization measurement 
is also not feasible with the present data sets. Precise v 2 and polarization 
measurements will need RHIC II luminosity. 
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Fig. 46. Left: The RHIC and SPS Raa data compared [175] with sequential char- 
monium suppression. The two dashed lines indicate no suppression (Raa/CNM = 
1) and complete suppression of if)' and Xc (Raa/CNM = 0.6). Right: The PHENIX 
midrapidity data compared with a model including a threshold density for J /if) 
suppression [149]. 

The sequential charmonium suppression model of Ref. [147] has been applied 
to the preliminary PHENIX data and to the SPS data. Normal nuclear absorp- 
tion is parameterized by an effective absorption cross section that accounts 
for all cold nuclear matter effects. The SPS p + A data give a larger effective 
J /if) absorption cross section than the RHIC d+Au data, as also implied by 
Ref. [114]. Three values of <r a b s are extracted from the d+Au data, one for 
each rapidity bin, similar to the ad hoc model applied to the PHENIX data 
in Figs. 39 and 40. These values are used to obtain the survival probability in 
A + A collisions. When N pax t is converted to energy density, e, and the survival 
probabilities for color screening and cold nuclear matter are included, the SPS 
and RHIC data were found to lie on a common suppression curve as a function 
of energy density. However the more precise, final PHENIX data, which show 
that the forward rapidity J /if) yield is significantly more suppressed than the 
midrapidity yield, are inconsistent with such a picture, as can be seen on the 
left-hand side of Fig. 46 [175]. The measured (pj) is also inconsistent with this 
scenario. 

The threshold energy density model applied to RHIC, Ref. [149], is compared 
to the midrapidity PHENIX data on the right-hand side of Fig. 46 for several 
values of the critical density. The model behavior is quite similar to that of the 
measured midrapidity Au+Au Raa for a critical density of 4 fm~ 2 . However, 
since it does not have any rapidity dependence, it cannot describe the much 
stronger suppression for Au+Au at forward/backward rapidity. It also cannot 



75 



explain the centrality dependence of the Cu+Cu suppression. 

The hydrodynamic model of Ref. [150] fits the critical temperatures to the 
midrapidity PHENIX Au+Au data but disagrees with R AA as a function of 
centrality. Since many effects are not included, it is difficult to draw conclu- 
sions. 

5.4 Proposed RHIC II quarkonia measurements 

Unlike other probes, quarkonia measurements are guided by predictions from 
lattice QCD calculations. Color screening modifies the linear rise of the QCD 
potential at large distances. The quarkonia spectral functions quantify the 
temperature dependence of the potential. Since quarkonia suppression is de- 
termined by the plasma temperature and the binding energy (equivalently the 
quarkonium size and the Debye screening length), measuring the sequential 
disappearance of these states acts as a QCD thermometer. 

Thus the importance of a comprehensive study of all experimentally acces- 
sible quarkonium states cannot be overstated. A systematic study of heavy 
quarkonium spectroscopy, with a complete determination of the suppression 
pattern of the quarkonium states, remains the most direct probe of decon- 
finement. It is also the signature that most closely resembles a thermometer 
of the hot initial state which, with future improved lattice calculations, can 
be directly compared to QCD. 

While J /ip physics is as compelling as it was in 1986 when first proposed 
by Matsui and Satz [27], the systematic study of all quarkonia states, and 
especially bottomonium, feasible at RHIC II, provides a more complete QGP 
probe than heretofore possible. 

Table 12 relates the main physics topics to the relevant probes and subsequent 
detector requirements. The ability of a program at RHIC II to make these 
measurements can be judged from the yields given in Tables 2,3, and 4. The 
measurements that are possible at RHIC without the luminosity upgrade are 
the J ftp rapidity and px distributions at full energy. The measurements that 
are newly possible at RHIC II are those of the excited charmonium states 
(ip f and Xc) and the bottomonium states (T(15 ( ), T(2S') and T(3S)). High 
Pt J/ip measurements, precise measurements of the J/ip v 2 and polarization, 
and excitation functions of heavy flavor measurements will be possible only 
at RHIC II. The possible statistical precision is illustrated in Fig. 47 which 
shows the statistical significance of the J/i/j nuclear modification factor Raa in 
Au+Au at high px, and Fig. 48, which presents the expected J/ip v 2 precision 
in Au+Au collisions at RHIC II [180]. It is evident that a comprehensive 
program to use quarkonium as a QCD thermometer to provide direct evidence 
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of deconfinement is possible only with RHIC II luminosity. 
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Fig. 47. The statistical significance after background subtraction of the J/tp 
Raa(pt) in PHENIX at central and forward rapidity from a 12 week Au+Au run 
with RHIC II luminosity (22 nb~ l delivered). Values above 5 GeV/c are a px ex- 
trapolation of current PHENIX data. 



The measurements needed to study the excited charmonium states, Xc and 
ip', have quite different problems. The ij}' measurement technique is the same 
as that for the J/ip, namely reconstruction of dilepton decays, but requires 
~ 100 times greater integrated luminosity for the same yield. In addition, the 
if)' measurement is more difficult because the background under the peak in 
the invariant mass spectrum is significant, increasing the integrated luminos- 
ity needed for precision measurements. The presence of the VTX detector in 
PHENIX will lead to significantly better mass resolution in the central arms 
by allowing a measurement of the full bend angle of the electrons - only the 
track angle after the magnetic field is measured at present. The FVTX de- 
tector is crucial to the if}' measurement in the muon arms because it both 
reduces the combinatorial background and improves the mass resolution, as 
shown in Fig. 49. A if>' measurement is certainly feasible at RHIC II. The \ c 
measurement can be done with the \ c —>■ J/if> 7 channel, where the J /if) is 
reconstructed from dilepton decays and the photon is detected in an electro- 
magnetic calorimeter. While the yields are larger than for the if)', the need 
to form the \ c invariant mass by combining each J/ip candidate with a large 
number of photons means that combinatorial backgrounds will be quite large 
in Au+Au collisions. Thus the Xc measurement will be difficult in central 
heavy-ion collisions. 

PHENIX has very recently carried out detailed simulations of the performance 
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Physics Motivation 


Probes 


Measurements 


Requirements 


Baseline measurements 


JM V, Xc T(15), T(25), 
and T(3S') decays to dilep- 
tons 


Rapidity and p T spectra in 
p + A and pp as a function 

Of y/S NN 


High luminosity and accep- 
tance for sufficient statistics, 
especially for the T family. 
Good mass resolution to re- 
solve ip and T states. 


Deconfinement and initial 
temperature 


J/iP, V, Xc T(15), T(25), 
and T(3S) decays to dilep- 
tons 


A + A suppression patterns as 
a function of ^s^n and A 


High luminosity, accep- 
tance and mass resolution 
for quarkonium, and trig- 
gers that work in Au+Au 
collisions. 


Thermalization and trans- 
port 




J/ip t>2 as function of ^snn 
and A 


High luminosity for good 
statistics in short runs for 
^s NN and A scans. 



Table 12. The main physics goals of the RHIC II quarkonium program with corresponding probes, studies, and requirements. 



PHENIX, Au+Au RHIC II, 12 weeks 

cm 0.25 

> 

~5 0.2 



0.15 
0.1 
0.05 


0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

J/y p T (GeV/c) 

Fig. 48. The vertical bars show the estimated absolute precision of the J ftp V2 
measurement in PHENIX at central and forward rapidity in a 12 week Au+Au run 
with RHIC II luminosity (22 nb~ x delivered) [180]. The precision is compared with 
several different models [17,48]. 

of the PHENIX muon arms in combination with the Nose Cone Calorimeter for 
Xc —> J r /i> 7 measurements. Figure 50 shows the expected precision for the Raa 
of the Xc using the suppression pattern from a simple model in which the Xc 
disappears at 1.16T C . Note that the maximum N pait bin ends at 30% centrality. 
Extending the measurement beyond 30% to the most central collisions is still 
being worked on. 

As is the case for the J/if>, the bottomonium states are studied through their 
dilepton decays. The bottomonium measurements require very large integrated 
luminosity and good invariant mass resolution. PHENIX expects to be able to 
resolve the T(1S), T(2S) and T(3S) states, see Fig. 51. Because of its larger 
acceptance, STAR will have substantially greater T yields than PHENIX. To 
clearly resolve the three S-wave states, STAR needs to reduce the amount of 
material between the interaction point and the inner field cage of the TPC 
because energy-loss due to bremsstrahlung deteriorates the mass resolution. 
Figure 52 shows the T mass spectrum without the material from the inner 
tracking system. Whether the resolution with the new inner tracking system 
(HFT) installed is sufficient to resolve all 3 states depends on the final design 
parameters of the new detectors. 

Although the yields are small relative to the J/ip, bottomonium measure- 
ments are quite clean. The states are massive (~ 10 GeV/c 2 ) so that their 
decay leptons have relatively large momenta and are thus easily distinguished 
from background leptons. The combinatorial background is small and multiple 
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Fig. 49. The J ftp and ip' invariant mass spectrum in the PHENIX muon arms 
without (left) and with (right) the improvement in mass resolution from the FVTX 
detector. 
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Fig. 50. The expected precision of the Raa for the Xc from PHENIX using 
X c — > J ftp 7 decays measured in the muon arms and the Nose Cone Calorimeter. 
The yields shown are for a delivered Au+Au luminosity of 30 rib . 



scattering is of less concern. While the interpretation of charmonium suppres- 
sion is made more difficult by the rather large cross section for nucleon and 
comover absorption, the situation for bottomonium is considerably better. Ab- 
sorption of directly produced bottomonium by hadronic comovers was shown 
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to be negligible [181]. 
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Fig. 51. The dielectron mass spectrum for the T family from a PHENIX central 
arm simulation, showing the expected improvement in T mass resolution provided 
by the initial direction measurement in the SVTX barrel [54]. The number of events 
shown correspond to a delivered luminosity of 30 nb~ l for Au+Au. 



The ^s^n dependence of produced J ftp's relative to the number of cc pairs, 
depicted in Fig. 29 [140], is striking. Measurement of the excitation function for 
this ratio over 30 < ^/snn < 200 GeV could help to disentangle suppression 
from enhancement. Such measurements, however, are extremely demanding 
statistically since both heavy quarks and quarkonia will need to be measured 
with good statistics over a wide energy range. 

A measurement of the quarkonium nuclear modification factor at high can 
provide a unique experimental probe of energy loss and color diffusion [182] . At 
relatively large transverse momentum, suppression due to color screening and 
coalescence are predicted to be negligible. Instead, the quarkonium state is a 
hard probe that interacts with the medium. In particular, any color octet can 
suffer energy loss. The relative abundance of charmonium resonances can pro- 
vide an experimental handle on such phenomena as each resonance may have a 
different octet contribution. We must exercise caution, however, as competing 
charmonium production models exist. In parallel with nucleus-nucleus studies, 
it is therefore important to investigate and compare production mechanisms 
in pp and p + A interactions, at both central and forward rapidities [183-185]. 

In addition to the baseline quarkonium measurements in pp and p + A colli- 
sions listed in Table 12, other measurements are required as model input for 
description of the A + A measurements. The most prominent of these are listed 
in Table 13. 



81 



160 
140 
I 120 

100 

to 



r STAR 

Simulation of 
Y Resolution 



80 



.a 

CD 



2 60 
* 40 



20 



8.5 



Sum 

Y(1S) 

Y(2S) 

Y(3S) 

Resolution without material from 
Inner Silicon Tracker System 




9.5 10 
Invariant Mass (GeV/c 2 ) 



10.5 



11 



Fig. 52. The T family dielectron mass spectrum simulated by STAR. The three 
states can be clearly separated in the data given sufficient statistics. 



The importance of measuring the underlying charm distributions as input to 
models of J/ip coalescence is obvious, as is the importance of understanding 
cold nuclear matter effects on quarkonium production. 



It is crucial for the interpretation of the A + A quarkonia yields to understand 
the feed down contributions from the Xc states, see Fig. 31. The best feed down 
measurement will be made in 500 GeV pp collisions because the increased 
luminosity and increased charmonium production cross sections lead to ~ 10 
times larger charmonium yields than in 200 GeV pp collisions. Since the Xc 
contribution to the J/ip yield will not change significantly between 200 and 
500 GeV, the increased yield at 500 GeV will provide a definitive baseline 
measurement of Xc feed down in pp collisions. 



Recently, quarkonium polarization measurements were suggested as signatures 
of QGP formation [186]. The quarkonium yields at RHIC II will be large 
enough to permit a J/ip polarization measurement at low px by both PHENIX 
and STAR. 
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Topic 


Measurements 


Requirements 


Cold nuclear effects 


In pp and p+A collisions: 

• -Tin y 7-1 ann 11 n p- 

pendence of quarkonia 
production 

• A dependence 


Large rapidity accep- 
tance, including forward 
coverage. 


Suppression vs coales- 
cence 


In pp, p + A and A + A 
collisions: 

• Charm da/dpTdy 

• J/y> t> 2 

• Raa(pt) 


High resolution vertex 
detectors (charm). 


Feed down 


Xc, at least in pp and p + 


Photon detection over 
wide rapidity range. High 
rates, good energy and 
momentum resolution to 
enhance Xc signal to 
background. 


Production mechanism 


Xc, polarization at least 
in pp and p + A 


Large acceptance for 
cos 6* measurement. 



Table 13 



Baseline measurements (beyond Au+Au) required in order to address the main 
physics questions. 

6 Relationship to the LHC program 

The major differences between quarkonium studies at RHIC II and at the LHC 
will be the temperature and lifetime of the medium, the relative production 
cross sections, luminosities and run times. 

The initial temperature in ^snn = 5.5 TeV central Pb+Pb collisions at the 
LHC is expected to be ~ 4 T c , while it is ~ 2 T c in y / J/v r /v= 200 GeV central 
Au+Au collisions at RHIC [187]. The QGP lifetime at the LHC is expected 
to be two to three times longer than at RHIC [187]. 

Heavy flavor production cross sections are much larger at the LHC. The open 
charm and bottom production cross sections are ~ 15 and ~ 100 times higher 
respectively [15]. The charmonium and bottomonium cross sections are ~ 13 
and ~ 55 times higher respectively than at RHIC [50] . The higher LHC open 
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heavy flavor cross sections increase the number of cc and bb pairs produced 
in central A + A collisions. There are ~ 10 and ~ 0.05 pairs, respectively, 
in Au+Au collisions at RHIC, while there should be ~ 115 and ~ 5 pairs in 
central Pb+Pb collisions at the LHC [15]. 

The RHIC II Au+Au average luminosity is projected to be 10 times larger than 
the LHC Pb+Pb luminosity (5 x 10 27 cm" 2 s _1 relative to 5 x 10 26 cm" 2 s" 1 ). 
The yearly heavy-ion runs at RHIC II are also expected to be considerably 
longer than at the LHC. Taking the polarized pp program at RHIC II into 
account, the heavy-ion program is expected to get ~ 12 week physics runs on 
average per year while the LHC heavy ion program will be allocated a one 
month physics run per year. Thus the annual integrated luminosity at RHIC 
II is expected to be about 30 times larger for heavy ions than at LHC. 

The larger heavy flavor cross sections at the LHC are approximately balanced 
by the increased luminosity and running times at RHIC II, making the heavy 
flavor yields per year similar. Thus the types of measurements that can be 
made at the two facilities will also be similar as well as of similar quality 
(see Tables 2, 3 and 5). However, there will be important differences in the 
physics environments prevailing at the two facilities which will make the two 
programs complementary. 

The higher initial energy density at the LHC means that the QGP will be 
created at a significantly higher temperature with a correspondingly strong 
potential for new physics effects at the LHC. In addition, the factor of ten 
increase in cc pairs and the factor of 100 increase in bb pairs per central 
collision at the LHC will have a major impact on the interpretation of heavy 
flavor measurements. We will discuss some of those differences here. 

Lattice calculations suggest that the J/ip may remain bound at the highest 
RHIC temperatures, while the excited charmonium states are predicted to be 
unbound. At the LHC, all of the charmonium states should be unbound at the 
highest temperatures, implying that almost all charmonium production in cen- 
tral Pb+Pb collisions at the LHC will be due to coalescence of cc pairs. Thus 
the prompt charmonium yields at the LHC should reflect only the coalescence 
mechanism with no contribution from the primordial J/ip production except 
in very peripheral collisions. The measurements at RHIC and the LHC will 
thus provide very different windows on charmonium suppression in the QGP 
that will help resolve the ambiguities in interpreting data due to the balance 
between destruction and coalescence formation of charmonium at RHIC. 

Because of its higher binding energy, the characteristics of bottomonium pro- 
duction at the LHC should be similar to those of charmonium at RHIC. The 
bottomonium states are shown in Fig. 53. The T(1S) may remain bound at 
the highest temperatures at the LHC while the other bottomonium states will 
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be dissociated. Given ~ 5 bb pairs in central Pb+Pb collisions (relative to 
~ 10 cc pairs at RHIC), the T yield at the LHC is predicted [188] to reflect 
a balance between dissociation and coalescence reminiscent of the RHIC J/ip 
production models. However, at RHIC, the bottomonium dissociation rates 
will be significantly different. While the T(1S) is predicted to be bound, the 
T(2S) may also remain bound. Only the T(3S) is likely to dissociate at RHIC. 
Also, since the bb pair yield at RHIC is ~ 0.05 per central Au+Au collision, 
no significant bottomonium production by coalescence is expected. Thus the 
bottomonium yields at RHIC II should reflect only QGP suppression. Mea- 
surements at RHIC II and the LHC will thus provide very different windows 
on bottomonium suppression in the QGP that will help to resolve the ambi- 
guities in interpretation due to the balance of bottomonium destruction and 
coalescence at the LHC. 

T (11020) 
T (10860) 



T(4S) 

BB threshold 



V3S) 

hadrons 



V2S) 




T(1S) 

V 1S > 

jPC = o-+ 1— 1 + - 0++ 1++ 2++ 

Fig. 53. Bottomonium mass levels and spin states. The common feed down channels 
are indicated. 

The open heavy flavor programs at RHIC II and the LHC will consist of 
similar measurements with similar goals. They will study energy loss, ther- 
malization and flow of heavy quarks in systems with very different energy 
densities, interaction cross sections and lifetimes. However, not all challenges 
in the measurements are similar. At ^snn = 200 GeV, bottom decays to 
leptons begin to dominate the single electron spectrum at px ~ 4 GeV/c. As 
the collision energy increases, the lepton spectra from B and D decays move 
closer together rather than further apart. Thus, the large increase in the bb 
cross section relative to cc does not make single leptons from B and D decays 
easier to separate. Preliminary calculations show that the B — > e decay does 
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become larger than that of D — > e but only for pt > 10 GeV/c. The two lep- 
ton sources differ by less than a factor of two to pt ~ 50 GeV/c in the range 
\y\ < 1. Separating single leptons from charm and bottom decays will require 
statistical separation using differences in the displaced vertex distributions at 
all pr at the LHC. Thus interpretation of single lepton data from heavy flavor 
decays will be more difficult at the LHC. 

ALICE can reconstruct D° decays from px ~ to pr ~ 25 GeV/c [51]. 
Like STAR, ALICE will be unable to trigger on D°'s and will have to obtain 
these events from the minimum bias sample. Thus the longer running times 
at RHIC are an advantage since more minimum bias data can be taken (see 
Tables 3 and 5). While it is not yet clear what CMS and ATLAS will do 
to reconstruct charm, they should be able to do b jets well, similar to the 
Tevatron measurements. As at RHIC, B mesons can be measured cleanly at 
the LHC through their decays to J/tp, although triggering on low p T J/ip's is 
difficult at the LHC. 

It has also been suggested that the BB contribution to the dimuon continuum, 
the dominant contribution above the T mass, can be used to measure energy 
loss [52] . That channel would be fairly clean at the LHC but more difficult at 
RHIC. 



7 Conclusions 



We have shown that, so far, the RHIC heavy flavor physics program has been 
very rich and stimulating, with many provocative and challenging results. To 
fully realize the potential of this compelling program, however, both detector 
upgrades and a luminosity upgrade are mandatory. Detector upgrades will 
improve reconstruction of charm hadron decays into hadronic channels and 
allow detection of B — > J/ipX decays using secondary vertex measurements. 
Upgrades will also make Xc detection possible and, in the case of STAR, lead to 
significant quarkonium yields. However only increased luminosities will allow 
high statistics measurements of all of these yields as well as increase the pr 
reach of J/ip and heavy flavor R A a and v 2 - 

We have also shown that the RHIC II and LHC heavy flavor physics pro- 
grams are complementary. Both are required for a complete understanding of 
heavy flavor production as a function of energy and temperature. We have 
also demonstrated that, despite lower heavy flavor cross sections at RHIC, 
the longer running times and higher luminosity of RHIC II make the recorded 
yields similar at the two facilities. 



86 



Acknowledgements 



We thank F. Karsch for the lattice contribution. The work of R. Vogt was 
supported in part by the U. S. Department of Energy under Contract Nos. 
DE-AC52-07NA27344 (LLNL) and DE-AC02-05CH11231 (LBNL) and by the 
National Science Foundation Grant PHY-0555660. The work of T. Ullrich 
was supported in part by the U. S. Department of Energy under Contract No. 
DE-AC02-98CH10886. The work of A. D. Frawley was supported by National 
Science Foundation grant PHY-04-56463. 



References 



[I] I. Arsene et al. [BRAHMS Collaboration], Nucl. Phys. A 757 (2005) 1 
[arXiv:nucl-ex/0410020] . 

[2] B. B. Back et al. [PHOBOS Collaboration], Nucl. Phys. A 757 (2005) 28 
[arXiv:nucl-ex/0410022]. 

[3] J. Adams et al. [STAR Collaboration], Nucl. Phys. A 757 (2005) 102 
[arXiv:nucl-ex/0501009] . 

[4] K. Adcox et al. [PHENIX Collaboration], Nucl. Phys. A 757 (2005) 184 
[arXiv:nucl-ex/0410003]. 

[5] M. Gyulassy and L. McLerran, Nucl. Phys. A 750 (2005) 30 [arXiv:nucl- 
th/0405013]. 

[6] M. Cacciari, P. Nason and R. Vogt, Phys. Rev. Lett. 95 (2005) 122001 
[arXiv:hep-ph/0502203]. 

[7] R. Vogt, Eur. Phys. J. Special Topics 155 (2008) 213 [arXiv: 0709. 2531 [hep- 
ph]]. 

[8] R. Vogt, Phys. Rev. C 71 (2005) 054902 [arXiv:hep-ph/0411378]. 

[9] R. Vogt, Acta Phys. Hung. New Ser. Heavy Ion Phys. 25 (2006) 97 [arXiv:nucl- 
th/0507027]. 

[10] M. Djordjevic et al, Phys. Lett. B 632 (2006) 81 [arXiv:nucl-th/0507019]. 

[II] M. Djordjevic, M. Gyulassy, and S. Wicks, Phys. Rev. Lett. 94 (2005) 112301 
[arXiv:hep-ph/0410372]. 

[12] N. Armesto, C. A. Salgado, and U. A. Wiedemann, Phys. Rev. D 69 (2004) 
114003 [arXiv:hep-ph/0312106]. 

[13] Yu. L. Dokshitzer and D. E. Kharzeev, Phys. Lett. B 519 (2001) 199 
[arXiv:hep-ph/0106202]. 



87 



[14] Z.-W. Lin, R. Vogt and X.-N. Wang, Phys. Rev. C 57 (1998) 899 [arXiv:nucl- 
th/9705006]; Z.-W. Lin and R. Vogt, Nucl. Phys. B 544 (1999) 339 [arXiv:hep- 
ph/9808214]. 

[15] R. Vogt [Hard Probe Collaboration], Int. J. Mod. Phys. E 12 (2003) 211 
[arXiv:hep-ph/01 11271]. 

[16] B. Svetitsky and A. Uziel, Phys. Rev. D 55 (1997) 2616 [arXiv:hep- 
ph/9606284]. 

[17] V. Greco, C. M. Ko and R. Rapp, Phys. Lett. B 595 (2004) 202 [arXiv:nucl- 
th/0312100]. 

[18] Z.-W. Lin and D. Molnar, Phys. Rev. C 68 (2003) 044901 [arXiv:nucl- 
th/0304045]. 

[19] B. I. Abelev et al. [STAR Collaboration], Phys. Rev. Lett. 98 (2007) 192301 
[arXiv:nucl-ex/0607012] . 

[20] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 98 (2007) 172301 
[arXiv:nucl-ex/0611018]. 

[21] S. Baumgart [STAR Collaboration], arXiv:0709.4223 [nucl-ex]. 

[22] R. L. Thews and M. L. Mangano, Phys. Rev. C 73 (2006) 014904 [arXiv:nucl- 
th/0505055]. 

[23] R. L. Thews, M. Schroedter, and J. Rafelski, Phys. Rev. C 63 (2001) 054905 
[arXiv:hep-ph/0007323]. 

[24] A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, Phys. Lett. B 
571 (2003) 36 [arXiv:nucl-th/0303036]. 

[25] A. P. Kostyuk, M. I. Gorenstein, H. Stocker and W. Greiner, Phys. Rev. C 68 
(2003) 041902 [arXiv:hep-ph/0305277]. 

[26] M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 410 (1997) 327. 

[27] T. Matsui and H. Satz, Phys. Lett. B 178 (1986) 416. 

[28] F. Karsch, M. T. Mehr, and H. Satz, Z. Phys. C 37 (1988) 617. 

[29] S. Datta, F. Karsch, P. Petreczky and I. Wetzorke, J. Phys. G 30 (2004) S1347 
[arXiv:hep-lat/0403017]. 

[30] D. Kharzeev and H Satz, Phys. Lett. B 334 (1994) 155 [arXiv:hep- 
ph/9405414]. 

[31] P. Braun-Munzinger and J. Stachel, Phys. Lett. B 490 (2000) 196 [arXiv:nucl- 
th/0007059]. 

[32] P. Braun-Munzinger and J. Stachel, Nucl. Phys. A 690 (2001) 119 [arXiv:nucl- 
th/0012064]. 



88 



[33] L. Grandchamp and R. Rapp, Phys. Lett. B 523 (2001) 60 [arXiv:hep- 
ph/0103124]. 

[34] P. Braun-Munzinger and K. Redlich, Eur. Phys. J. C 16 (2000) 519 [arXiv:hep- 
ph/0001008]. 

[35] C. M. Ko, B. Zhang, X.-N. Wang and X. F. Zhang, Phys. Lett. B 444 (1998) 
237 [arXiv:nucl-th/9808032]. 

[36] J. F. Gunion and R. Vogt, Nucl. Phys. B 492 (1997) 301 [arXiv:hep- 
ph/9610420]. 

[37] N. Armesto, M. Cacciari, A. Dainese, C. A. Salgado and U. A. Wiedemann, 
Nucl. Phys. A 774 (2006) 589 [arXiv:hep-ph/0510284]; Phys. Lett. B 637 
(2006) 362 [arXiv:hep-ph/0511257]. 

[38] A. Mischke [STAR Collaboration], in proceedings of Quark Matter 2008. 

[39] Y. Morino [PHENIX Collaboration], in proceedings of Quark Matter 2008. 

[40] A. Adare et al. [PHENIX Collaboration], arXiv:0802.0050 [hep-ex]. 

[41] D. Hornback [PHENIX Collaboration], in proceedings of Quark Matter 2008. 

[42] A. Adare et al. [PHENIX Collaboration] Phys. Rev. Lett. 98 (2007) 232301 
[arXiv:nucl-ex/0611020]. 

[43] A. Adare et al. [PHENIX Collaboration] arXiv:0801.0220vl [nucl-ex]. 

[44] A. Adare et al. [PHENIX Collaboration] Phys. Rev. Lett. 98, 232002 (2007) 
arXiv:hep-ex/0611020vl. 

[45] A. Capella and E. G. Ferreiro, Eur. Phys. J. C 42 (2005) 419 [arXiv:hep- 
ph/0505032]. 

[46] L. Grandchamp, R. Rapp and G. E. Brown, Phys. Rev. Lett. 92 (2004) 212301 
[arXiv:hep-ph/0306077]. 

[47] A. Adare et al. [PHENIX Collaboration], arXiv:0711.3917vl [nucl-ex]. 

[48] L. Yan, P. Zhuang, and N. Xu, arXiv:nucl-th/0608010. 

[49] X. Zhao and R. Rapp, arXiv:0712.2407 [hep-ph]. 

[50] M Bedjidian et al, arXiv:hep-ph/0311048. 

[51] A. Dainese, R. Vogt, M. Bondila, K. J. Eskola and V. J. Kolhinen, J. Phys. G 
30 (2004) 1787 [arXiv:hep-ph/0403098]. 

[52] I. P. Lokhtin and A. M. Snigirev, Eur. Phys. J. C 21 (2001) 155 [arXiv:hep- 
ph/0105244]. 

[53] O. Kodolova, CMS Internal Note 2003-002 (2003). 

[54] PHENIX Silicon Vertex Detector proposal, PHENIX Collaboration (2004) 
BNL Report BNL-72204R-2004. 



89 



[55] PHENIX Forward Vertex detector proposal, PHENIX Collaboration (2007) 

http://www.phenix.bnl.gov/phenix/WWW/docs/upgrades/fvtx/PHENIX_FVTX- 
proposal_ll-l-07.pdf 

[56] PHENIX Nose Cone Calorimeter proposal, PHENIX Collaboration (2006) 

http:/ /www. phenix.bnl.gov/phenix/WWW/docs/upgrades/ncc/PHENIX_NCC_march2006.pdf. 

[57] J. Lajoie [PHENIX Collaboration], AIP Conf. Proc. 915 (2007) 490 
[arXiv:physics/0612018]. 

[58] K. H. Ackermann et al, Nucl. Instrum. Meth. A 499 (2003) 624. 

[59] W. J. Llope, Nucl. Instrum. Meth. B 241 (2005) 306. 

[60] G. Van Nieuwenhuizen [STAR Collaboration] (2007) private communication. 

[61] R. V. Gavai, D. Kharzeev, H. Satz, G. A. Schuler, K. Sridhar and R. Vogt, 
Int. J. Mod. Phys. A 10 (1995) 3043 [arXiv:hep-ph/9502270]. 

[62] B. Alessandro et al. [ALICE Collaboration], J. Phys. G: Nucl. Part. Phys. 32 
(2006) 1295. 

[63] David G. d'Enterria, (Ed.) et al. J. Phys. G 34 (2007) 2307. 

[64] L. Rosselet (ATLAS Collaboration) CERN-ATL-PHYS-CONF-2007-011, HEP 
conference 2007. 

[65] M. Arneodo, Phys. Rept. 240 (1994) 301. 

[66] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 94 (2005) 062301 
[arXiv:nucl-ex/0407006] . 

[67] R. Vogt, M. Cacciari and P. Nason, Nucl. Phys. A 774 (2006) 661. 

[68] M. Cacciari, M. Greco and P. Nason, JHEP 9805 (1998) 007 [arXiv:hep- 
ph/9803400]; M. Cacciari, S. Frixione and P. Nason, JHEP 0103 (2001) 006 
[arXiv:hep-ph/0102134]. 

[69] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B 303 (1988) 607; 
Nucl. Phys. B 327 (1989) 49 [Erratum: Nucl. Phys. B 335 (1990) 260]. 

[70] W. Beenakker et al, Nucl. Phys. B 351 507. 

[71] M. Cacciari and M. Greco, Nucl. Phys. B 421 (1994) 530 [arXiv:hep- 
ph/9311260]. 

[72] M. Cacciari, private communication. 

[73] M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 122003 [arXiv:hep- 
ph/0204025]. 

[74] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B 592 
(2004) 1. 

[75] A. Adare et al. [PHENIX Collaboration] Phys. Rev. Lett. 97 (2006) 252002 
[arXiv:hep-ex/0609010]. 



90 



[76] J. D. Bjorken, FERMILAB-PUB-82-059-THY (1982). 

[77] M. H. Thoma and M. Gyulassy, Nucl. Phys. B 351 (1991) 491 

[78] M. H. Thoma, Phys. Lett. B 273 (1991) 128. 

[79] S. Mrowczynski, Phys. Lett. B 269 (1991) 383. 

[80] E. Braaten and M. H. Thoma, Phys. Rev. D 44 (1991) 2625. 

[81] B. Svetitsky, Phys. Rev. D 37 (1988) 2484. 

[82] Y. Koike and T. Matsui, Phys. Rev. D 45 (1991) 3237. 

[83] M. G. Mustafa, D. Pal, D. K. Srivastava and M. Thoma, Phys. Lett. B 428 
(1998) 234 [arXiv:nucl-th/9711059]. 

[84] E. Shuryak, Phys. Rev. C 55 (1997) 961 [arXiv:nucl-th/9605011]. 

[85] M. Djordjevic and U. Heinz, arXiv:0802.1230 [nucl-th]. 

[86] S. Wicks, W. Horowitz, M. Djordjevic and M. Gyulassy, Nucl. Phys. A 784 
(2007) 426 [arXiv:nucl-th/0512076]. 

[87] A. Adil and I. Vitev, Phys. Lett. B 649 (2007) 139 [arXiv:hep-ph/0611109]. 

[88] G. D. Moore and D. Teaney, Phys. Rev. C 71 (2005) 064904 [arXiv:hep- 
ph/0412346]. 

[89] H. van Hees and R. Rapp, Phys. Rev. C 71 (2005) 034907 [arXivmucl- 
th/0412015]. 

[90] H. van Hees, V. Greco and R. Rapp, Phys. Rev. C 73 (2006) 034913 
[arXiv:nucl-th/0508055]. 

[91] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. D 76 (2007) 092002 
[arXiv:hep-ex/0609032]. 

[92] S. Kelly [PHENIX Collaboration], J. Phys. G 30 (2004) S1189. 

[93] C. Zhong [STAR Collaboration], J. Phys. G 34 (2007) S741 [arXivmucl- 
ex/0702014]. 

[94] A. A. P. Suaide, J. Phys. G 34 (2007) S369 [arXiv:nucl-ex/0702035]. 

[95] R. A. Lacey et al, Phys. Rev. Lett. 98 (2007) 092301 [arXiv:nucl-ex/0609025]. 

[96] S. Gavin and M. Abdel-Aziz, Phys. Rev. Lett. 97 (2006) 162302 [arXivmucl- 
th/0606061]. 

[97] B. V. Jacak and M. P. McCumber, J. Phys. G 34 (2007) S543 [arXivmucl- 
ex/0703032]. 

[98] V. D. Barger, W. Y. Keung, and R. J. Phillips, Phys. Lett. B 91 (1980) 253. 
[99] V. D. Barger, W. Y. Keung, and R. J. Phillips, Z. Phys. C 6 (1980) 169. 



91 



[100] G. A. Schuler and R. Vogt, Phys. Lett. B 387 (1996) 181 [arXiv:hep- 
ph/9606410]. 

[101] M. L. Mangano, P. Nason and G. Ridolfi, Nucl. Phys. B 405 (1993) 507. 
[102] R. Vogt, arXiv:hep-ph/0203151. 

[103] R. Vogt, Acta Phys. Hung. New Ser. Heavy Ion Phys. 18 (2003) 11 [arXiv:hep- 
ph/0205330]. 

[104] A. D. Martin et al, Eur. Phys. J. C 4 (1998) 463 [arXiv:hep-ph/9803445]. 

[105] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 12 (2000) 375 
[arXiv:hep-ph /9903282] . 

[106] M. Gluck, E. Reya and A. Vogt, Eur. Phys. J. C 5 (1998) 461 [arXiv:hep- 
ph/9806404]. 

[107] S. Digal, P. Petreczky and H. Satz, Phys. Rev. D 64 (2001) 094015 [arXiv:hep- 
ph/0106017]. 

[108] T. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 84 (2000) 2094 
[arXiv:hep-ex/9910025]. 

[109] S. R. Klein and R. Vogt, Phys. Rev. Lett. 91 (2003) 142301 [arXiv:nucl- 
th/0305046]. 

[110] D. Kharzeev and H. Satz, Phys. Lett. 366 (1996) 316 [arXiv:hep-ph/96xxxx]. 

[Ill] R. Vogt, Nucl. Phys. A 700 (2002) 539 [arXiv:hep-ph/0107045]. 

[112] K. J. Eskola, V. J. Kolhinen and P. V. Ruuskanen, Nucl. Phys. B 535 (1998) 
351 [arXiv:hep-ph/9802350]; K. J. Eskola, V. J. Kolhinen and C. A. Salgado, 
Eur. Phys. J. C 9 (1999) 61 [arXiv:hep-ph/9807297]. 

[113] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 96 (2006) 012304 
[arXiv:nucl-ex/0507032]. 

[114] M. Leitch, C. Lourenco and R. Vogt, J. Phys. G 34 (2007) S759. 

[115] L. McLerran and B. Svetitsky, Phys. Lett. B 98 (1981) 195; J. Kuti, J. Polonyi 
and K. Szlachanyi, Phys. Lett. B 98 (1981) 199. 

[116] O. Kaczmarek and F. Zantow, Eur. Phys. J. C 43 (2005) 63 [arXiv:hep- 
lat/0502011]. 

[117] S. Digal, P. Petreczky and H. Satz, Phys. Lett. B 514 (2001) 57 [arXiv:hep- 
ph/0105234]. 

[118] C. Y. Wong, Phys. Rev. C 72 (2005) 034906 [arXiv:hep-ph/0408020]. 

[119] W. M. Alberico, A. Beraudo, A. De Pace and A. Molinari, Phys. Rev. D 72 
(2005) 114011 [arXiv:hep-ph/0507084]. 

[120] E. V. Shuryak and I. Zahed, Phys. Rev. D 70 (2004) 054507 [arXiv:hep- 
ph/0403127]. 



92 



[121] D. Blaschke, O. Kaczmarek, E. Laermann and V. Yudichev, Eur. Phys. J. C 
43 (2005) 81 [arXiv:hep-ph/0505053]. 



[122 
[123 
[124 

[125 
[126 

[127; 
[128; 

[129 
[130 
[131 
[132 
[133 
[134 
[135 
[136 

[137; 



[138 
[139 



H. Satz, J. Phys. G 32 (2005) R25 [arXiv:hep-ph/0512217]. 

F. Karsch, Eur. Phys. J. C 43 (2005) 35 [arXiv:hep-lat/0502014]. 

Y. Nakahara, M. Asakawa, and T. Hatsuda, Phys. Rev. D 60 (1999) 091503 
[arXiv:hep-lat /9905034] . 

N. Brambilla et al. [Quarkonium Working Group], arXiv:hep-ph/0412158. 
F. Karsch and R. Petronzio, Z. Phys. C 37 (1988) 627. 

M. Asakawa and T. Hatsuda, Phys. Rev. Lett. 92 (2004) 012001 [arXiv:hep- 
lat/0308034]. 

S. Datta, F. Karsch, S. Wissel, P. Petreczky and I. Wetzorke, arXiv:hep- 
lat/0409147. 

K. L. Haglin and C. Gale, Phys. Rev. C 63 (2001) 065201 [arXiv:nucl- 
th/0010017]. 

A. Mocsy and P. Petreczky, Phys. Rev. Lett. 99 (2007) 211602 
[arXiv:0706.2183 [hep-ph]]. 

T. Umeda, R. Katayama, O. Miyamura and H. Matsufuru, Int. J. Mod. Phys. 
A 16 (2001) 2215 [arXiv:hep-lat/0011085]. 

H. Iida, T. Doi, N. Ishii and H. Suganuma, PoS LAT2005 (2006) 184 
[arXiv:hep-lat/0509129]. 

M. I. Gorenstein, A. P. Kostyuk, H. Stoecker and W. Greiner, Phys. Lett. B 
509 (2001) 277 [arXiv:hep-ph/0010148]. 

M. I. Gorenstein, A. P. Kostyuk, H. Stoecker and W. Greiner, J. Phys. G 
27(2001) L47 [arXiv:hep-ph/0012015]. 

M. I. Gorenstein, A. P. Kostyuk, L. McLerran, H. Stoecker and W. Greiner, 
arXiv:hep-ph/0012292. 

R. L. Thews, Nucl. Phys. A 702 (2002) 341 [arXiv:hep-ph/0111015]. 

K. Martins, D. Blaschke, and E. Quack, Phys. Rev. C 51 (1995) 2723 
[arXiv:hep-ph/9806027]; S. G. Matinyan and B. Miiller, Phys. Rev. C 58 
(1998) 2994 [arXiv:nucl-th/9806027]; K. Haglin, Phys. Rev. C 61 (2000) 
031902R [arXiv:nucl-th/9907034]; Z.-W. Lin and C. M. Ko, Phys. Rev. C 62 
(2000) 034903 [arXiv:nucl-th/9912046]. 

K. Geiger and J. R. Ellis, Phys. Rev. D 52 (1995) 1500 [arXiv:hep-ph/9503349]; 
Phys. Rev. D 54 (1996) 1967 [arXiv:hep-ph/9511321]. 

R. L. Thews, Eur. Phys. J. A 29 (2006) 15 [arXiv:hep-ph/0511292]. 



93 



[140] L. Grandchamp and R. Rapp, Nucl. Phys. A 715 (2003) 545 [arXiv:hep- 
ph/0209141]. 

[141] S. Gavin, M. Gyulassy and A. Jackson, Phys. Lett. B 207 (1988) 257. 

[142] R. Vogt, M. Prakash, P. Koch and T.H. Hansson, Phys. Lett. B 207 (1988) 
263. 

[143] J. Ftacnik, P. Lichard and J. Pisut, Phys. Lett. B 207 (1988) 194; J. Ftacnik, 
P. Lichard, N. Pisutova and J. Pisut, Z. Phys. C 42 (1989) 132. 

[144] A. Capella et al, arXiv:0712.4331vl [hep-ph]. 

[145] N. Armesto, A. Capella and E. G. Ferreiro, Phys. Rev. C 59 (1999) 395 
[arXiv:hep-ph/9807258]. 

[146] A. Capella and E. Ferriero, arXiv:hep-ph/0610313. 

[147] F. Karsch, D. Kharzeev and H. Satz, Phys. Lett. B 637 (2006) 75 [arXiv:hep- 
ph/0512239]. 

[148] J.-P. Blaizot and J.-Y. Ollitrault, Phys. Rev. Lett. 77 (1996) 1703 [arXiv:hep- 
ph/9606289]. 

[149] A. K. Chaudhuri, Phys. Rev. C 75 (2007) 044902 [arXiv:nucl-th/0610031]. 

[150] T. Gunji, H. Hamagaki, T. Hatsuda and T. Hirano, Phys. Rev. C 76 (2007) 
051901 [arXiv:hep-ph/0703061]. 

[151] H. Liu, K. Rajagopal, and U. A. Wiedemann, Phys. Rev. Lett. 98 (2007) 
182301 [arXiv:hep-ph/0607062]. 

[152] A. L. S. Angelis et al. [HELIOS-3 Collaboration], Eur. Phys. J. C 5 (1998) 63. 

[153] M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 466 (1999) 408. 

[154] M. C. Abreu et al. [NA50 Collaboration], Phys. Lett. B 477 (2000) 28. 

[155] A. Baldit et al. [NA60 Collaboration], CERN-SPSC-2000-010 (2000). 

[156] H. Kolanoski [HERA-B Collaboration], J. Phys. G 31 (2005) S799. 

[157] I. Abt et al. [HERA-B Collaboration], Phys. Lett. B 561 (2003) 61 [arXiv:hep- 
ex/0211033]. 

[158] G. Borges [NA50 Collaboration], J. Phys. G 30 (2004) S1351. 
[159] B. Alessandro et al. [NA50 Collaboration], arXiv:nucl-ex/0612012. 
[160] C. Lourenco, R. Vogt and H. Woehri, in preparation. 

[161] M. Leitch et al. [FNAL E866/NuSea collaboration], Phys. Rev. Lett. 84 (2000) 
3256 [arXiv:nucl-ex/9909007]. 

[162] B. Alessandro et al. [NA50 Collaboration], Eur. Phys. J. C 49 (2007) 559 
[arXiv:nucl-ex/0612013] . 



94 



[163] B. Alessandro et al. [NA50 Collaboration], Eur. Phys. J. C 39 (2005) 335 
[arXiv:hep-ex/0412036]. 

[164] M. Nardi and H. Satz, Phys. Lett. B 442 (1998) 14 [arXiv:hep-ph/9805247]. 

[165] D. Kharzeev, C. Lourenco, M. Nardi and H. Satz, Z. Phys. C 74 (1997) 307 
[arXiv:hep-ph/9612217]. 

[166] R. Arnaldi et al. [NA60 Collaboration], arXiv:0706.4361 [nucl-ex]. 

[167] J.-P. Blaizot, M. Dinh and J.-Y. Ollitrault, Phys. Rev. Lett. 85 (2000) 4012 
[arXiv:nucl-th/0007020]. 

[168] R. Shahoyan [NA60 Collaboration], Eur. Phys. J. C 43 (2005) 209. 

[169] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 92 (2004) 051802 
[arXiv:hep-ex/0307019]. 

[170] Z. Tang [STAR Collaboration], in proceedings of Quark Matter 2008. 

[171] H. Buesching [PHENIX Collaboration], Nucl. Phys. A 774 (2006) 103; 
M. Leitch, J. Phys. G 34 (2007) S453 [arXiv:nucl-ex/0701021]. 

[172] P. Djawotho et al. [STAR Collaboration], J. Phys. G 34 (2007) S947 
[arXiv:nucl-ex/0701075]. 

[173] D. deFlorian and R. Sassot, Phys. Rev. D 69 074028 (2004) 

[174] B. Alessandro et al, Euro. Phys. J. C 48 329 (2006). 

[175] T. Gunji [PHENIX Collaboration], J. Phys. G 34 (2007) S749 [arXivmucl- 
ex/0703004]. 

[176] R. Rapp, Eur. Phys. J. C 43 (2005) 91 [arXiv:hep-ph/0502208]. 

[177] A. Andronic, P. Braun-Munzinger, K. Redlich and J. Stachel, Nucl. Phys. A 
789 (2007) 334 [arXiv:nucl-th/0611023]. 

[178] W. Cassing, E. L. Bratkovskaya, and S. Juchem, Nucl. Phys. A 674 (2000) 
249 [arXiv:nucl-th/0001024]. 

[179] X. Zhao and R. Rapp, arXiv:0712.2407v2 [hep-ph]. 

[180] M. Wysocki (2007), private communication. 

[181] Z.-W. Lin and C. M. Ko, Phys. Lett. B 503 (2001) 104 [arXivrnucl- 
th/0007027]. 

[182] R. Baier, D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50 (2000) 
37 [arXiv:hep-ph/0002198]. 

[183] R. Vogt, Phys. Rev. C 61 (2000) 035203 [arXiv:hep-ph/9907317]. 

[184] S. Gavin and J. Milana, Phys. Rev. Lett. 68 (1992) 1834. 

[185] M. B. Johnson et al, Phys. Rev. C 65 (2002) 025203 [arXiv:hep-ph/0105195]. 



95 



[186] B. L. Ioffe and D. E. Kharzeev, Phys. Rev. C 68 (2003) 061902 [arXiv:hep- 
ph/0306176]. 

[187] I. Vitev, J. Phys. G 30 (2004) S791 [arXiv:hep-ph/0403089]. 

[188] L. Grandchamp, S. Lumpkins, D. Sun, H. van Hees and R. Rapp, Phys. Rev. C 
73 (2006) 064906 [arXiv:hep-ph/0507314]. 



96 



